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INTRODUCTION 

' The promise of filamentary composite materials,, whose de- 

velopment may be considered as entering its second generation, 
continues to generate intense interest and applications activ- 
' ity. Such interest and activity are well-founded, since they 

are based on the possibility of using relatively brittle mate- f' 

I rials with high modulus, high strength, but low density in 

composites with good durability and high tolerance to damage 
and which, when they do fail, do so in a non-catastrophic 
manner. Fiber reinforced composite materials of this kind 
offer substantially improved performance and potentially 
lower costs for aerospace hardware. 

Much progress has been achieved since the initial devel- 
' opments in the mid 1960 *s. Rather limited applications to 

primary aircraft structure have been made, however, mainly in 
a material-substitution mode on military aircraft, except for 
a few experiments currently underway on large passenger air- 
planes in commercial operation and a few military develop- 
ments which have not seen service use. 

To fulfill the promise of composite materials completely 
requires a strong technology base. NASA and AFOSR recognize 
the present state of the art to be such that to fully exploit 
composites in sophisticated aerospace structures, the tech- 
nology base must be improved. This, in turn, calls for ex- 
panding fundamental knowledge and the means by which it can 
be successfully applied in design and manufacture. 

PRECEPING PAGE BLANK NOT FILMED 
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Aa technology of compoaito materials and structures moves 
toward fuller adoption into aerospace struaturos, some of the 
problems of an earlier era are being solved, others which 
seemed important are being put into porspoctivo as relatively 
minor, and still others unanticipated or put aside are emerg- 
ing as of high priority. The purpose of the RPI program as 
funded by NASA and AFOSR has been to develop critical advanced 
technology in the areas of physical properties, structural 
concepts and analysis, manufacturing, reliability and life 
prediction. 

Our approach to accomplishing those goals is through an 
interdisciplinary program, unusual in at least two important 
aspects for a university. First, the nature of the research 
is comprehensive - from fiber and matrix constituent proper- 
ties research, through the integration of constituents into 
composite materials and their characterization, the behavior 
of composites as they are used in generic structural compo- 
nents, their non-destructive and proof testing, to the logical 
conclusion of such activities; namely .research into the com- 
posite structure's long term integrity under conditions perti- 
nent to service use. Inherent in the RPI program is the moti- 
vation which basic research into the structural aspects pro- 
vides for research at the materials level, and vice versa. 

Second, interactions among faculty contributing to pro- 
gram objectives - which is a group wider than that supported 
under the project - is on a day to day basis, regardless of 


organizational lines. Program management is largely at the 
working level, and administrative, scientific and technical 
decisions are made, for the most part, independent of con- 
siderations normally associated with academic departments. 
Involvement of this kind includes - depending on the flow of 
the research - faculty, staff and students from chemistry, 
chemical engineering, civil engineering, materials engineer- 
ing, and the department of mechanical engineering, aeronau- 
tical engineering and mechanics. 

Both of these characteristics of the NASA/APOSR program 
of research in composite materials and structures foster the 
kinds of fundamental advances which are triggered by insights 
into aspects beyond the narrow confines of an individual dis- 
cipline. This is a program characteristic often sought in 
many fields at a university, but seldom achieved. 

Overall program emphasis, is on basic, long-term re- 
search in the following categories: (a) constituent mater- 

ials, (b) composite materials, (c) generic structural ele- 
ments, (d) processing science technology and (e) maintaining 
long-term structural integrity. Progress in the program will 
be reported in the following pages under these headings. 

Those computer software developments are also undertaken which 
both support Rensselaer projects in composite materials and 
structures research in the areas listed above and which also 
represent research with the potential of v/idely useful re- 
sults in their own right. 
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In short, the NASA/AFOSR Composites Aircraft Program is 
a multi-faceted program planned and managed so that scien- 
tists and engineer, 3 in a number of pertinent disciplines will 
Interact to achieve its goals. Research in the basic compo- 
sition, characteristics and processing science of composite 
materials and their constituents iis balanced against the 

mechanics, conceptual design, fabrication and testing of gen- 
// 

eiric structural elements typical of aerospace vehicles so as 
to encourage the discovery of unusual solutions to present 
and future problems. In the following sections, more de- 
tailed descriptions of the progress achieved in the various 
component parts of this comprehensive program are presented. 


PART II 

CONSTITUENT I4ATERIALS 


II-A THE SURFACE ENERGY OP ANHYDRIDE-CURED EPOXY RESINS 
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II-A THE SURFACE ENERGY OF ANHYDRIDE- CURED EPOXY RESINS 

Senior Investigator: R. J. Diefendorf 

1 « Introduction 

The comparatively low interlaminar shear strength of 
high modulus graphite composites compared to that of glass 
and boron fibers has resulted in an extensi.ve and intensive 
search for techniques to modify the fiber-matrix interface 
to optimize the interlaminar shear strength of these com- 
posites . 

Several ingenious routes have been taken to modify the 
fiber surfaces and to improve the bond between the fibers 
and epoxy matrices . Some of these approaches include oxi- 
dative treatments (Drzal*- ' •' , Donnet and Ehrburger^ •' ) 
which increase the shear strength of the composite but, 
generally, at the expense of composite, tensile strength. 
Interfacial condensation methods and electrolytic deposition 
processes which lead to the polymerization of monomers on 
the fiber surfaces have been developed to form a flexible 
interlink or tie coat. 

These various techniques have met with differing degrees 
of success. However, an effective and predictable method of 
modifying or designing the interlaminar shear strength of 
these composites is yet to be developed. The aims of this 
study are to describe the surface energy components of fibers 


Numbers in brackets in this section refer to the references 
which are listed on page 16 . 
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and resins and to further develop models which can predict 
these components from the architecture of the polymeric 
matrices . 

This research, conducted by graduate student C. E. 

Uzoh, is part of an effort to understand and catagorize the 
fundcimental principles which govern the design of an effec- 
tive and predictable interfacial bond between fiber and 
matrix. In this report, the semiempirical geometric mean 
model of Girifalco-Good-Fowkes-Kaeble-Young and Zisman was 
used to describe the role of anhydride curing agents in de- 
termining the surface energies of cured epoxy resins. 

2. Status 

During the last period, the surface energy components 
of various epoxy resins, curing agents and several non- 
epoxide polymers were characterized. 

3. Progress During Report Period 
a. Surface Energy Components 

V7etting studies of various liquids on the surfaces of 
resins indicate that the curing agent is the critical factor 
that dictates the dispersive and polar fractions of the sur- 
face energies of these epoxies . Epoxies studied include 
nadic methyl anhydride (NMA) and dodecenyl succinic anhy- 
dride (DDSA) -cured bisphenol A, novolac, cycloaliphatic and 
tetraglycidyl derivatives of epoxy with benzyldimethylamine 
(BDMA) as catalyst. 


For instance, the polar components for all DDSA-cured epoxies 
are consistently less than 5 erg/cm^, even though the total 
surface energy is about 30 erg/cm® (see Table II-A-1) . Also, 
DDSA-cured Epon 828, Epon 152, Aradi'te CY-179 and MY-720 have 
about the same value of and y , irrespective of the anhy- 
dride/resin (A/R) ratio. Even though the NMA cures have 
higher y than DDSA cures, the surface energy compo- 

nents have the same insensitivity to A/R ratio. This insen- 
sitivity of surface energy components to A/R ratio is due to 
the splitting-off of the excess epoxy components during post- 
curing at elevated temperature. 

Comparative studies on surfaces of amine-cured epoxies 
and other non-epoxides indicate that the observed difference 
in NMA and DDSA-cured samples is due to the architecture of 
the curing agents. The difference is caused by the flexi- 
bility induced by the long aliphatic chain and the steric 
hindrance of methylene groups in DDSA, compared to the more 
compact structure of NMA. This explanation is supported by 
the microhardness of the polished surfaces of these epoxies 
(see Table II-A-2) . 

b. Critical Surface Energy 

The critical surface energy y of a low energy solid 

o 

was operationally defined by Zisman as: 

Cos 6 = 1 + b(y^ - y^^) 

where 6 is the contact angle, y^^ the surface energy of the 
contacting liquid and b is a constant. Homologous series of 
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TABLE II-A-1 

SURFACE ENERGY COMPONENTS OF ANHYDRIDE -CURED EPOXIES 


Substrates 

Dispersive 

Component 

Polar 

Component 

Total Surface 
Energy* 

** 

DDSA Cure 





CY - 179 
(95-100 phr) 

29.2 ± 0.6 

2.9 ± 0.7 

32.1 ± 

1.4 

MY - 720 
(68-102 phr) 

26.2 ± 1.2 

3.5 ± 1.4 

29.8 ± 

2.6 

Epon - 828 
(130-140 phr) 

30.1 ± 0.6 

2.5 ± 0.7 

32.6 ± 

1.4 

Epon - 152 
(106-144 phr) 

27.1 + 1.2 

2.3 ± 1.5 

29.4 ± 

2.7 

NMA Cure 





CY - 179 
(62-74 phr) 

34.3 ± 1.2 

10.2 ± 1.4 

44.4 ± 

2.6 

MY - 720 
(52-^6 phr) 

35.5 ± 0.7 

9.2 ± 0.8 

44.6 ± 

1.5 

Epon - 152 
(77-89 phr) 

34.3 ± 0.8 

7.8 ± 1.0 

42.1 ± 

1.8 

DTA/Epon - 828 
(12 phr, 24°C) 

36.4 + 0.7 

8,5 ± 0.8 

42.9 ± 

1.5 

* . ? 

All values in erg/ cm 

DDSA (dodecenyl succinic anhydride) 

*** 




NMA (nadic methyl anhydride) 


ff 
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TABLE II"A"2 

MICROHARDNESS OF ANHYDRIDE-CURED EPOXIES 

Microhardnesa (khn) 

i( ifii 

Resins NMA Cure DDSA Cure 

Aradite CY-179 
Aradite MY-720 
Epon - 152 
Epon - 828 

NMA (nadic methyl anhydride) 

^ A 

DDSA (dodecenyl succinic anhydride) 


22.0 ± 2.7 
24.7 ± 2.3 
25.6 ± 2.9 


12.6 ± 1.0 
14.4 ± 1.2 


12.4 ± 0.4 


liquids are known to obey this approximate modal, however, 

where non-homologous series of liquids are used, the value 

of y for a given solid is poorly defined because of scatter, 

f7] 

J. L. Gordon ■■ observed that y of polymers tends to 
increase with the solubility parameter, fi, and he related Y 

c 

to the surface energy of a solid, y, as; 

Y = Yc/<l>^ 

where 0 refers to the interaction between the polymer and 
the class of liquids used in the experiment. 

r g y| 

L. H. Lee^ ' / by taking the entropy term into consid- 

eration in the modified he at-of -mixture equation, related y 
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to the cohesive energy density, 6^, and the molar volume, 
in the form of: 


5 2 = 16,81- 


Y /V 
'c' m 


1/3 


o.oc 


<!> = 


where (J) is a correction term for the deviation from the geo- 
metric mean model. Although y. is still a useful parameter 

c 

in characterizing the surface of a solid, its full meaning 
still remains a matter of conjecture. 

In this work, the critical surface energies, y^, of all 
the various epoxy and non-epoxide polymers studied, were al- 
ways observed to be somewhere between the dispersive com- 
ponent of the surface energy and the total surface energy 
predicted by their geometric mean model. 

Statistical analysis of the results indicated that, to a 
first approximation, y can be expressed as a linear combi- 

w 

nation of the dispersive and polar components of the surface 
energy (see Figure II-A-1) . This relationship can be con- 
veniently expressed as: 

Y^ = + ‘lY^ = 0.9759 

, N = 43 

Lim Yj, y 

y^ 0 


It can be argued that dispersive interactions will always be 
fully operable and the full contribution of the dispersive 
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component of the surface energy is required. For a random 
polymer solid with a constituent polar group, free rotation 
is severely restricted. Hence, polar-polar interaction will 
be considerably less than that of less restricted dipoles, 
as in liquid-liquid or liquid-gel interactions. In the 
liquid phase of a polymer, any reorientation by one molecule 
will invariably involve the cooperation of neighboring mole- 
cules or chain ends. Thus, the consequent reduction in 
polar interaction is an entropic effect. These screening 
effects and the inability of the molecules at the solid sur- 
face to reorient themselves to the desired configuration in 
response to the contacting liquid are best represented at 
low temperatures, when only half the value of the nondisper- 
sive contribution is used. 

The observed scatter in the experimental result is con- 
sistent with the standard deviations of the surface energy 
components of the liquids used in this study. From these re- 
sults, a fund£unental understanding of the interactions of 
fiber/resin interfaces is emerging. We hope to exploit this 
understanding during the next two report periods to determine 
if better interfaces can be designed. 

4 , Plans for Upcoming Period 

Future work will be devoted to a systematic study of 
the interfacial bond between various low energy surfaces, to 
graphite fiber surface modification and to techniques to meas- 
ure energy of fiber bundles. 
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6« Current Publications or ProBentations by 
Professor Diofondorf on this Subject 

"Physical, Thermal and Chemical Parameters Affecting the De- 
velopment of Pitch Mesophase", with S-H. Chen, S. L, Eilen- 
berg and D. M. Higgs. 

Published in Carbon '82, Society of the Chemical 
Industry, London, 1982, pp. 180-184. 

Presented at Carbon '82, London, UK, September 
23, 1982. 

"Discotic Liquid Crystals" 

Presented at CPRC Meeting, Bonbanne, France, 

September 17, 1982. 
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III--A ADVANCED ANALYSIS METHODS 

Senior Investigator: E. J. Brunelle 

1. Introduction 

The guest continues to find solutions for composite 
plates in terms of a minimum number of generalized variables. 
This approach continues to promise that a general parametric 
understanding of composite plates can be promulgated. 

2 . Status 

The principal results of the last period are summarized 
as follows: (a) The relations governing the location of 

double frequencies were developed for a simply-supported 
orthotropic plate; (b) The general solution, in terms of the 
similarity variable presented for the static de- 

formation, vibration and buckling of antisymmetric cross-ply 
plates with 52 boundary conditions and (c) The extremum solu- 
tions (D* = 0) for antisymmetric cross-ply plates were out- 
lined for arbitrary boundary conditions, and a sample solu- 
tion was presented, namely, the minimum uniaxial buckling 
coefficient envelope versus affine aspect ratio {k^ versus 
a^/b^) for antisymmetric cross-ply plates with all four 
edges clamped. 


3. Progress During Report Period 
Progress in the last reporting period has been made 
principally in the following two areas . 


th 
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First, the complete Karman-Rostovstev plate equations 
(including the in-plane and stress-strain laws) have been 
correctly restretched and several approximate solutions (as 
suggested by Donnell for the corresponding isotropic Karman 
equations) have been carried out as shown in Figures III-A- 
1 through 4 . 

Second, nevr material parameters R and X, defined as 

= (Ej/E2>/(Vj^2) ^ 

have been used to express the fundamental plate parameters 
D*, e and H*. These results are given by 

D* = i [1 + 2(1 - X"^)R] 

E = [1 + 2(1 - X”^)R]"* 

and using the values for fourteen materials collected from 
three sources (Tsai, Jones and Pipes) are plotted in Figure 
III-A-5 as three mini-figures. The range of the linear hori- 
zontal X scale is 0 £ X ^ 30 for all three mini-figures; the 
vertical one-cycle log scales for D*, e and 1/H* go from .1 
to 1.0 and the R parameter has values shown from .1 to 10. 
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4. Plans for Upcoming Period 
A general solution scheme will be attempted for cross- 
ply symmetric laminated plates , and an attempt will be made 
to extend the Karman-Rostovstev stretched equations for 
orthotropic plates to include cross-ply symmetric laminated 
plates. Also, extremum solution details will be scrutinized 
carefully in an effort to find similarity solutions for the 
orthotropic Mindlin plates (i.e. , shear deformable plates) . 

5. Current Publications or Presentations by 
Profeissor Brunelle on this Subject 

"Eigenvalue Similarity Rules for a Class of Rectangular 
Specially Orthotropic Laminated Plates" 

Presented at the Ninth U, S. National Congress of 
Applied Mechanics, Ithaca, NY, June 21-25, 1982. 

"A Similarity Variable for Generating the Frequency of Ini- 
tially Stressed Specially Orthotropic Laminates" 

To be presented at the AHS Composite Structures 
Specialist's Meeting, Philadelphia, PA, March 
23-25, 1983. 

"The Affine Equivalence of Local Stress Distributions in 
Damaged Composites and Batdorf's Electric Analogue" 

To be submitted to the AIAA Journal, December 1982. 
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III-B FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator; E, KrempX 

1. In t roduc tion 

The deformation and failure behavior of graphite-epoxy 
tubes uvider biaxial (tension, torsion) loading is being in- 
vestigated with the assistance of Ph.D. candidate T.-M. Niu. 
The aim of this research is to provide basic understanding 
and design information on the biaxial response of graphite- 
epoxy composites. 


2 . Status 

Using graphite-epoxy [±45] thin-walled tubes, the 
static axial and torsional strength were determined together 
with the elastic moduli. Fatigue tests under completely re- 
versed, load-controlled axial loading showed comparatively 
poor axial fatigue strength; furthermore, fatigue life was 

affected by frequency, A paper on this subject has appeared 

rn* 

in the Journal of Composite Materials ■' . 

3. Progress During Report Period 
During the report period, the previously generated data 
were analyzed and interpreted. The analyses are contained in 
the forthcoming Doctor of Engineering thesis of Mr. Niu's. 

Of particular interest was the observation of time-de- 
pendent effects in axial and combined axial-torsion cyclic 
loading. 

*Numbers in brackets in this section refer to the references 
which are listed on page 35. 
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a. Tlmo-Depondence of Deformation 

The deformation behavior of the [±45]^ thin-walled 
graphite-epoxy tubes was found to be characterized by a fiber- 
controlled torsional and a matrix-controlled axial behavior. 

It was shown that significant time-dependent deformation 
ensued at stress levels exceeding twenty-five percent of the 
axial ultimate strength while, in torsion, time dependent ef- 
fects were negligible. 

The studies of time-dependent deformation during biaxial 
loading can be obtained from in-phase multistep creep and re- 
laxation tests. These tests were performed in an MTS servo- 
controlled testing machine controlled by an MTS 463 Data/ 
Control Processor. 

After stabilizing the properties of the specimen by 
applying two completely reversed triangular loading cycles, 
multistep creep or relaxation loadings were then performed 
in the next two cycles under load or strain control in in- 
pliase loading; i.e., with maximum torsion applied at the 
same instant as maximum axial load. The amplitude of the 
oscillatory load or displacement was held constant at 0, 50 
and 100 percent of the maximum value to be used in the tests. 
The frequency was chosen to be 0.05 Hz. The hold time for 
each step was 9 cycles or 180 seconds. 

In the case of in-phase loading, the stress component 
angle, G, in the o-x plane was chosen to be 131°. The negative 
twist then coincides with the tensile stretch, thus eliminating 


the torsional buckling. The test results show that local buck 
ling, in fact, did not occur (see the T-y diagram in Figure 
XIl-B-1) . In both hysteresis loops, the most pronounced 
creep strains wore observed at the stress level extremes of 
cycles 3 and 4. It was observed that the creep strains at 
zero stress level are not the smallest encountered. Due to 
fiber dominance of deformation, the torsional creep strain at 
eaph step is much less than the corresponding axial strain. 

By switching the control mode from stress to strain in 
the same running program, multistep in~phase relaxation can 
be studied. The strain component angle, 0', in e-y plane was 
chosen to be 162®. The hysteresis loops au a frequency of 
0.05 Hz for torsion and axial deformation are shown in Figure 
III-B-2. The maximum relaxed stresses. At and Aa, were found 
at the extrema of the strain levels. Again, the relaxed 
shear stress at each step is less than the corresponding re- 
laxed axial stress. 

In-phase, load-controlled loading applied along 6 =49® 
constitutes the most severe loading for the occurrence of 
local buckling. In this case the negative twist and the com- 
pression occur simultaneously. The hysteresis loops for such 
a loading are shown in Figure III-B-3. The creep effects in 
the tensile part developed sim. larly to those shown in Figure 
III-B-1, and the torsional creep strains for positive torque 
are generally small. When, at cycle 3, creep hold times were 
introduced, local torsional buckling was observed at the 


I 





Figure III-B-l. Hysteresis Loops during In-phase Kultlstep Creep, Three-Minute Periods 
of Creep, Ffeq. *0.05 Hz, 0 •» 131“ in the a-x Plane (Tube 88). 
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minimum load level - as evi.denced by the increase of shear 
strain (see points 2 and 3 in Figure lll~B-3) . On subse- 
quent cycles the torsional creep strain repeated itself. This 
example shows that time-dependent effects Can be found even 
under fiber controlled deformation, when combined loading is 

I 

imposed. 


4. Plans for Upcoming Period 
The properties of tubes with a new lay-up [0/±45] of 

5 

graphite-epoxy material will be determined. In this case, 
axial and torsional deformation are fiber dominated. The 
fabrication of several batches of this new specimen is in 
progress. Subsequently, the static and fatigue properties 
of this tube will be explored. 

5 . References 

1. Krempl, E. and T.-M. Niu, Journal of Composite Materials, 
Vol 16, 1982, pp. 172-187. 

6 . Current Publications or Presentations by 
PJ^ofQssor Krempl on this Subject 

"Biaxial In-Phase and Out-of-Phase Behavior of Graphite- 
Epoxy Tubes" 

To be presented at the ASTM Conference on Biaxial/ 
Multiaxial Fatigue, C-i-n Francisco, California, 

December 15-17, 1982. 


37 


III-C VISCOELASTIC CHARACTERIZATION OF IN SITU RESINS AND 
NEAT RESINS 

Senior Investigator; S. S. Sternstein 

1 . Introduction 

This project is concerned with those properties of high 
performance composites which are strongly dependent on the 
physical properties of the matrix resin. To date, the re- 
search has involved the precise viscoelastic characterization 
of epoxy neat resins, interlaminar failure of composites and 
the inhomogeneous swelling of and the effects of moisture on 
composites . 


2. Status 

The effects of moisture uptake in neat resins and com- 
posites on the viscoelastic characteristics (loss factor, 
dispersion anomolies and location of the glass transition of 
the matrix) of various samples has been investigated. The 
composite samples have been tested using the centro-symmetric 
deformation (CSD) geometry described previously. Samples were 
cut from laminates having a (0/90) stacking sequence and 
prepared from prepreg material s plied by the Fiberite Cor- 
poration (HYE1048A1E) . Neat resin samples were tested in 
uniaxial tension using thin rectangular strips. 

3. Progress During Report Period 
A set of composite samples exposed for various periods 
of time to a boiling water environment were tested (prior to 
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any drying history) and the in-j.’jhase (storage) stiffness (M*) 
results are shown in Figure III-C-1 versus weight gain mois- 
ture) . The values of M' are normalized by the value of M' 
for the dry samiple and also normalized by the cube of the 
thickness for the wet sample (as tested) . It is clear that 
at about 2.5% moisture pickup, the slope of the curve marked- 
ly increases , suggesting a change in the structure-moisture 
interaction mechanism. 

The results of a similar set of experiments on neat 
resins is shown in Figure IlI-C-2, v;here the slope change is 
also obvious, but occurs at about 7% moisture. Since the 
composite samples are about 35% resin, this suggests that 
the slope change in the composite is due to the excessive 
softening of the matrix and is, therefore, matrix controlled. 
The variation of loss factor (M'*/M* the so-called tan 6) 
with weight gain for the neat resin samples is given in Fig- 
ure III-C-3. The ordinate represents the increase (or incre- 
ment) in loss factor which is obtained in a wet sample rela- 
tive to a dry sample . 

During the previous report period, we reported that 
large permanent increases in loss factor and decreases in 
stiffness occur when a wet composite sample is dried and then 
tested. This effect (hysteresis) does not occur, or is mini- 
mal, in the neat resin samples. Thus, in Figure III-C-3 a 
dried sample would lie on the curve shown and not exhibit a 
large residual loss factor. For the reader's convenience. 
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Figure III-C-4 is a repeat of our previous result on loss 
factor hysteresis for coinposites. 

The stoichiometric relationship between the stiffness- 
moisture relationship of the composite CFigure III-C-1) and 
the neat resin (Figure III-C-2) , the absence of drying hys- 
teresis on the loss factor of the neat resin (Figure III-C-3) 
and the presence of hysteresis in the composite (Figure III- 
C-4) suggest that (a) the matrix stiffness is a significant 
factor in the out-of-plane stiffness of the composite; (b) 
matrix stiffness is very nonlinear with respect to moisture 
content and a significant mechanism change occurs at about 
7% moisture content (resin basis - 948) and (c) drying hys- 
teresis in the composite represents a failure mechanism either 
between the fibers and matrix and/or delamination. 


4 . Plwns for Upcominc Period 


It xs planned to conclude this study of moisture uptake 
and drying effects on neat resins and composites and to sup- 
plement the mechanical data with suitable microscopic evi- 
dence of the failure patterns. 
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III-D Nl!I*lE?aCAL INVESTIGATION OF MOISTURE EFFECTS 

Senior Investigator; M. S, Shephard 

1 . Introduction 

This analysis builds on the theoretical developments of 
Professor S. Sternstein (reported in progress reports dated 
December ’79 and June ’80) in which the fully-coupled, non- 
linear thermomechanical equations of inhomogeneous swelling 
of composites in the presence of moisture and temperature 
are analyzed. Here the one-dimensional case of a single 
fiber in an infinite matrix has been extended to two-dimen- 
sional, multiple fiber cases with far field stresses. 

2. Status 

As discussed in earlier progress reports, moisture ef- 
fects were introduced into the problem through a nonlinear 
constitutive relation, the principle of virtual displacements 
was employed to develop a nonlinear matrix equation using 
displacement-based finite elements and these nonlinear equa- 
tions were solved using quasi-Newton methods. Subsequently, 
the constitutive equation was built into the finite element 
analysis. The results of a single fiber comparison analysis 
and two fiber analysis were also given. 

3. Progress During Report Period 

During the last reporting period, graduate student Frida 
Lumban-Tobing has continued to run addition problems and to 

PRECEDING PAGE BLANK NOT FILMED 
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collect data on both the recults obtained and the computer 
solution times required using different quasi-Newton updates. 
The additional cases run include a four-fiber square pack 
arrangement and a seven-fiber hexagonal pack arrangement. 

Both cases considered one radius fiber spacing under 100% 
humidity conditions, with and v/ithout far field stresses. 

a. Exa m ple of Problems 

a.l) Square Packed Four Fiber Case . Figure III-D-1 
shows the basic model considered for the case of four fibers 
in a square pack at a one radius spacing while Figure III-D-2 
shows the finite element mesh of six-noded quadratic isopara- 
metric elements for the quarter of the plate that was con- 
sidered. The cases analyzed for this configuration included 
zero far field stress and various combinations of uniform 
non-zero far field stress, ali at 100% humidity. 

There is a strong interaction between the fiber and 
matrix for the zero load case, as indicated in Figure III- 
D-3, which shows a set of maximum principal stress contours 
in the area directly around the fiber. The peak stress val- 
ues in the matrix for this case are approximately 25% higher 
than for the unloaded single fiber case presented previously. 

Three different loading cases were considered for this 
problem under 100% humidity. The first was a uniform load 
in one direction (PX) , uniform loads of the same sign in 
both directions (PX = PY) and opposite uniform loads in both 
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directions (PX « -PY) . Figure III-D-4 summarizes the re- 
sults for these load cases. This figure plots the maximum 
value of the first principal stress obtained at any point in 
the matrix as a function of load for the various load cases. 
The trends are similar to those obtained for the two fiber 
case with the highest stresses occurring for the case of 
PX ^ “PY. 


a. 2) Hexagonal Packed Seven Fiber Case . Figure III-D-5 
shows the basic model considered for this case using a mini- 
mum fiber spacing of one radius. The mesh of six-noded 
quadratic isoparametric triangles used to analyze this case 
is shown in Figure III-D-6. The cases analyzed for this 
configuration also include zero far field stress and various 
combinations of non-zero far field stress at 100% humidity. 

Again there is a strong interaction between the fiber 
and matrix for the zero load case, as indicated in Figure 
III-D-7, which shows a set of maximum principal stress con- 
tours in the area around the fiber. The peak stress value 
in the matrix for this case is 10% higher than the one fiber 
case . 

Four different 100% humidity loading cases were consid- 
ered for this problem. They are uniform x-direction load 
(PX) , uniform y-direction load (PY) , uniform loads of the 
same sign in the x- and y-direction (PX = PY) and uniform 
loads of opposite signs in the two directions (PX = -PY) . 
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Figure III-D-8 summarizes the results for these four load 
cases. This figure plots the maximum value of the first 
principal stress obtained at any point in the matrix as a 
function of load level . These trends are similar to those 
obtained for the two fiber and four fiber cases. 

b. Comparison of Solution Techniques 

The previous progress report discussed the various 
quasi-Newton update methods used in the program. During 
this reporting period, a study was done to determine the 
relative efficiency of the Davidon, Broyden and BFGS up- 
dates. Table III-D-1 gives a comparison of the number of 
operations required in the calculation of the increment step 
at the k iteration and the addition storage required past 
the initial triangulation for the three update methods and 
a standard Newton method. As can be seen from this table, 
the Newton method requires far more operations per step than 
any of the quasi-Newton updates. 

Consideration of the convergence characteristics. Table 
III-D-2, indicates that the quasi-Newton methods require much 
less computer time than a standard Newton method, with the 
Davidon update requiring less than one-half the time. A 
comparison of the three updates indicates that the Davidon 
update required the least number of iterations, fewer line 
searches and the least computer time. Thus, it is concluded 
that the Davison update is the most efficient for this prob- 


lem. 
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OPERATION COUNTS AND ADDITIONAL STORAGE 
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mk(k 3) 


mk 


2k + 4 + 



k 

k(2) 



Calculation of increment step at k 


th 


iteration 


m = Number of equations 


Additional 

Storage 


k(m 4* 1) 

k(2m) 

k(2m) 


n - Bandwidth 


ORIGINAL PAGE IS) 
OF POOR QUALITY 


56 


TABLE III--D--2 
CONVERGENCE COMPARISON 


Problem 


Method CPU* 


No. of 
Iterations 


No, of 
Searches 


Hexagonal Pack 

Davidon 

17:21 

9 

1 

P = 0 

Broyden 

20:26 

9 

5 

H = 100% 

BFGS 

21:32 

11 

4 

start ='>•'>= 

Newton 

38:40 

8 

6 

Square Pack 

Davidon 

10:55 

10 

3 

= 7.5 

H = 100% 

~o ” ^ p = 0 

Broyden 

14:34 

14 

5 

BFGS 

14:48 

16 

3 

Two Fibers 

Davidon 

6:51 

10 

6 

(distance lOR apart) 

Broyden 

15:07 

19 

19 

P = 0 

BFGS 

9:13 

11 

11 


H = 50% 


start 


0.024 


Times in minutes and seconds on 


a PRIME 


500 


computer 
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4. Plans for Upcoming Period 
During the next reporting period, Frida Lumban-Tobing 
will be finishing the final analysis studies and concluding 
her doctoral research on this topic. 


5. Current Publications or Presentations by 
Professor Shephard on this Subject 

"Finite Element Analysis of Moisture Effects in Graphite- 
Epoxy Composites", with F. E. Lumban-Tobing and S. S. Stern- 
stein. 

To be published in Advances and Trends in Structural 
and Solid Mechanics, A. K, Noor and J. M. Housner, 

Eds., October 1982 and Computers and Structures, 

1983. 

To be presented at the Symposium on Advances and 
Trends in Structural and Solid Mechanics , Washing- 
ton, D. C. by F. E. Lumban-Tobing, October 4, 1982. 
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III-E NUMERICAL INVESTIGATION OF THE MICROMECHANICS OF 
COMPOSITE FRACTURE 

Senior Investigator; M. S. Shephard 
I, Introduction 

To understand the mechanisms of failure in composites 
it is necessary to develop insight into the micromechanical 
behavior, including interactions between matrix and fibers 
as the load is increased from zero to that corresponding to 
failure. Investigations of these phenomena, either experi~ 
mental or numerical, are difficult. The purpose of this 
project, being carried out by graduate student Nabil Yehia, 
is to develop a nonlinear finite element analysis capability 
for performing numerical investigations of the micromechan- 
ical failure of composites. 

2 . Status 

The overall structure of a general purpose finite ele- 
ment code for fracture problems has been designed and is 
currently being implemented and tested. The major design 
considerations for the program are that it can accomodate 
various types of fracture problems and house a variety of 
algorithmic approaches to fracture and debonding analysis . 
The linear portion of the program has been coded and tested. 
The initial fracture criteria have been selected and are 
currently being implemented. 
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3* Progress During Report Poriod 
During the last reporting period, effort has been con- 
centrated on implementing and testing the initial computa- 
tional models for fracture mechanics computations. The areas 
under current consideration are the introduction of the sin- 
gular elements, the determination of first cracking in an 
isotropic material, calculation of the crack tip stress in- 
tensity factors and the determination of the direction of 
crack propagation. 

g, singular Elements 

This is the simplest of the fracture criteria to imple- 
ment, since proper elastic crack tip singular elements can 
be obtained from the standard quadratic isoparametric ele- 

n 21 * 

ments by the proper positioning of the node points ' . 

Thus, a module has been placed in the program that will po- 
sition the nodes of all elements at the crack tip in the 
proper position. Tests have been run on both the singular 
quadratic isoparametric triangle and collapsed quadratic 
isoparametric quadrilateral. The triangular element checked 
out very well, while the collapsed quadrilateral demon- 
strated some minor problems in satisfying equilibrium. This 
problem is currently being investigated. 


Numbers in brackets in this section refer to the references 
which are listed on page 65. 
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b. Determination of First Cracking 

As mentioned in the previous progress report, the pro- 
gram is being designed to consider the possibility of a crack 
initiating under the action of the applied load. Thus far 
the maximum tangential stress criteron has been implemented 
and checked. This criterion predicts that a crack will start 
when the maximum tangential stress exceeds a specified value. 

c. Calculation of Stress Intensity Values 

The stress intensity factors and are being com- 
prted by the displacement method v/hich uses the local dis- 
placements along the crack surface and the two singular ele- 
ments on that surface. The major advantage of this approach 
is that its local nature makes it straightforrmrd to use in 
a variety of situations. 

To test this module, the example problem suggested in 

Reference [4] was analyzed. The geometry and dimensions used 

in the analysis are shown in Figure III-E-1. This problem 

is a pure mode I situation with an accepted formula for the 

r41 

stress intensity factor . Figure III-E-2 shows the geome- 
try analyzed and conditions applied. The numerical results 
obtained indicate that good results can be obtained for a 
reasonable number of elements. However, seemingly minor 
changes to the mesh layout can play an important role in de- 
termining the accuracy of the result. 

Figures III-E-3 and -4 show a typical element mesh used 
to analyze the problem. The error in obtained using this 
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Figure Geometry of Double Cracked Plate 


Edge Fixed Against 
Horizontal 
Motion . 


Xtg»- 


\ 


I 

mm 






Y Y Y Y Y 
Edge Fixed Against Verticle Motion 




Figure llI-E-2. 


Basic Idealization Used for the Finite Element 
Analysis of Double Cracked Plate 
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Figure III-E-3. Finite Element Mesh in Middle Section 

of Double Cracked Plate 
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Figure III-E-4. Close-up of Finite Element Mesh in 

Area for Double Cracked Plate 


1 

A 

< 

J 


63 


j 


ORlGSriAL PAPI- IG 
OF POOR o”ai;ty 


i 

< J 

i 

I 


I 


1 J 











basic Jjnosh configuration, but changing the grading of the 
elements along the crack, ranged from 18.6% to 0.19%. The 
reason for this variation has yet to be confirmed, but two 
factors are considered important. One of the factors is the 
well-known influence of ratio of the length of singular ele- 
ment edge to the length of the crack. The other, which is 
currently being considered, is the shape of the singular 
elements. These aspects will be investigated further, and 
the results obtained will be considered in the development 
of future meshes to insure accurate results. 

d. Determination of Direction of Crack Growth 

The purpose of this module is to determine the direc- 
tion of crack growth after it has been determined that the 
crack will propagate. The minimum strain energy factor cri- 
ter ion ^ ^ is being applied in this prediction. This lias 

been implemented and tested on the example problem discussed 
immediately above. The proper results were obtained, how- 
ever, further testing is required. 

4. Plans for Upcoming Period 
During the next reporting period, effort will concen- 
trate on continued implementing and testing of the various 
modules required for crack initiation and propagation prob- 
lems. In addition, some specific micromechanical analyses 
will be carried out on composites that can be properly repre 
sented with those features available in the program. 
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INTRODUCTION 

The redesign of flight-critical structural components 
currently in service as metal parts on existing aircraft was 
included in earlier stages of RPI's composite materials and 
structures program to assure relevance and realism to the 
project's direction and to contribute to the nation's pool 
of innovative structural design concepts. Two such compo- 
nents were selected in cooperation with major commercial 
airframe manufacturers, Boeing Commercial Airplane Company 
and Lockheed California Company, and with advice from NASA. 
The first was the 727 airplane's elevator actuator attachment 
rib, and the second was the wing-mounted engine drag strut 
for the L-1011 airplane. In previous progress reports, sev- 
eral areas of research in whigh results were needed for such 
redesigns - but were not then available - were identified. 
From among those, several which continue to represent gaps in 
needed technology are under investigation. Results to date 
will be reported, and probable directions of future work will 
be described. 
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IV~A MECHANICAL JOINTS IN COMPOSITES 

The 727 elevator actuator attachment rib redesign in 
composites at RPI, reported in earlier progress reports, led 
to separate research projects. One deals with heavily load- 
ed holes in composites such as are encountered at the actu- 
ator attachment; the second deals with lightly loaded mechan- 
ical joints, such as occur, for example, in a row of small 
bolts or rivets, typified by the more or less continuous 
joint between rib flange and the elevator skin. Both of 
these problems have required the development of new analysis 
and/or test methods . 

The conceptual redesign of the L-1011 drag strut has 
emphasized the importance and difficulty of designing com- 
pact lugs in composite materials. Two efforts pertinent to 
compact lug design have been pursued. First, unique struc- 
tural configurations have been investigated that may increase 
the strength of such lugs. Second, failure criteria appro- 
priate to the stress fields that exist in such lugs have 
been examined. 

All four of these projects are described in the para- 
graphs to follow. 
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1. Analysis of Heavily Loaded Mechanical Joints 

Senior Investigator: C. Muser 

a. Introduction 

Independent of the materials used, there often will be 
discontinuities in structural members, such as openings in 
panels and holes to accomodate mechanical joints. Although 
the fibers in advanced composites can be arranged in any 
imaginable manner (magnaweave, raagnaswire, stylusweave, 
etc.), the materials currently used for plates including 
those containing discontinuities are mainly orthotropic with 
respect to cartesian coordinates. The present investigation 
deals with cylindrically orthotropic materials, since they 
prordse to yield lower stresses around circular openings. 


b . Status 

A closed form solution for the stress distribution in a 
circular plate with a circular hole at its center, subjected 
to uniaxial traction, was found by N. J. Hoff for the case 
when its material is homogeneous and cylindrically ortho- 
tropic. The geometry of that plate can be seen in Figure 
IV-C-l-a of the 41st Semi-Annual Report. This analysis was 
extended to incorporate nonhomogeneous materials . The con- 
stitutive equations with compliances which are a function of 
the radius are: 


E . 
1 
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’'re'= ®66<1 + «66=''>''re 

where i,j = r, and ( ) indicates no sum. 

The exponent h describes the rate of change of material 
properties along the radius r. The constants Sfj and Q. . 
are a function of the size of the plate and the parameter 
h. The compliances S^j can be prescribed at both boundaries 
independently of the choice of the parameter h, so long as 
the convergence criteria are observed. 

Such material property variations allow the choice of 
structural characteristics which yield low stresses in the 
vicinity of discontinuities, but leave the capabilities of 
the plate as a structural element unaffected. The results 
were encouraging, and it was decided that the same analysis 
be used to treat a bolt-loaded circular plate. The geometry 
of such a plate is shown in Figure IV-C-l-b of the 41st 
Semi-Annual Report. The goal, in this case, is to develop 
an algorithm capable of handling any boundary conditions 
which are expressible in a Fourier series. To find the 
exact bondary conditions was considered beyond the scope of 
this phase of the research. 

c. Progress During Report Period 

The form of the solution selected incorporates stress 
functions, J, expressed as an infinite series. A typical 
stress function ^ is shown as follows; 
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There are four boundary conditions expanded into a series of 
n terms (n -»•«') . The corresponding constants are A^. For 
every n, the compatibility equation needs to be satisfied in 
its limit. The constants K™(n) emerge from this process. 

The homogeneous case yields the roots t^(n). 

A computer program was written to calculate the result- 
ing stresses. It needs as input the material properties at 
the boundary of the hole, at the outer boundary of the plate 
and the value of the parameter h. At this writing, the 
stress distribution in a plate subjected to uniaxial uniform 
traction (the empty hole case) can be calculated. The actu- 
al computation is very efficient, using about one dollar's 
worth of computer time. 

The program has been extended to handle the case of the 
bolt-loaded plate. However, more analytical work is needed 
to determine the roots of the first term in the series 
(n = 1) . 


d. Plans for Upcoming Period 
The analytical work needs to be completed, and the re- 
sults incorporated in the otherwise finished computer pro- 
gram. Once that has been done, stress distributions in cir- 
cular, bolt-loaded plates - the materials of which are cylin- 
drically orthotropic and vary radially - will be calculated 
numerically. 
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e. Current Publications or Presentations by 
Dr, Muser on this Subject 

"Stress Concentration Factors for Cy] indrically Orthotropic 
Plates", with N. J. Hoff. 

Published in the Journal of Composite Materials, 

Vol. 16, July 1982, pp. 313-317. 

"Stress Concentrations in Cylindrically Orthotropic Plates 
with Radial Variation of the Compliances", with N. J, Hoff. 

To be presented at the Fourth Conference on Compos- 
ite Materials (ICCM-IV) , Tokyo, Japan, October 
25-28, 1982. 
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2. Pin-Loaded Holes in Uniform Composite Plates 

Senior Investigator: U. G, Loewy 

a . Introduction 

In performing stress analysis of mechanical joints with 
clearances between fastener holes and loading pins (or bolts) , 
one encounters uncertain boundary conditions at points along 
the fastener hole edges. There are similarities in this sit- 
uation with the more usual elastic contact problems having 
uniform contact boundaries. 

The single fastener joint, the simplest case, has been 
studied as the first step, in the hope that this may be ex- 
tended to several other cases of interest. If the motions 
of hole edge points on the pin-hole contact boundary can be 
prescribed properly, then the only remaining task is to de- 
termine the amount of contact. It is clear that the amount 
of pin-hole contact is a function of applied load, joint 
geometry and the laminate properties of the plate. As dis- 
cussed in a previous report, this relationship is non-linear, 
and as a result, the whole problem is non-linear. The re- 
sults reported here are the work of doctoral candidate 
Wonsub Kim. 

b. Status 

A system of linear constraint equations has been intro- 
duced in a finite element method (FEM) analysis to describe 
the motions of hole edge points. The derivation of these 
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constraint equations is based on the fact that any hole edge 
point within the pin-hole contact region must move with the 
surface of the loading pin (assumed rigid and frictionless) 
and remain somewhere on the surface at equilibriiun for a 
given applied load. The region of contact may be represented 

by 26 (Figure IV-A-2-a) and is initially undetermined. 

c 

Three ways to determine 0 for any given load level 

o 

were established and explained in detail in the last report. 
They were (a) the "semi- empirical method”, in which isochro- 
matic fringes obtained in a photoelastic experiment were 
correlated with the maximum shear strain distributions cal- 
culated by FEM with an assumed value of 0 ; (b) the "compu- 
ter graphics terminal method", which seeks the 0 value 
which produces the smoothest deformed shape of the hole edge 
and (c) the "minimum energy potential of applied load (or 
minimum strain energy) method", which searches for the value 
of 0 that results in minimum energy potential for the given 
load application. All of these methods require iterative 
FEM analysis to find the correct 0^. 

Among these techniques , the "minimum energy potential 
of applied load method" is the most convenient, since its 
physical (or mathematical) interpretation is quite clear and 
it requires only minimum numbers of iterations, roughly four 
or five, as discussed in the previous report. 

Using these methods, analyses of a single fastener 
mechanical joint with pin-hole clearance (Figure IV-A-2-a) 
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were performed for two types of laminate^ and the results 
were presented and discussed in the previous report. Among 
the results, it was pointed out that the non-linearity of 
the stress behavior was evident, even for a very small 
amount of pin-hole clearance (0.0015 inch) . Further, a 
slight difference in laminate type, for example [±45®/0®]^g 
versus [±30°/0®] , resulted in quite distinguishable differ- 

ences in stress behavior. It was also observed that stress 
increases, in individual layers with differing fiber orien- 
tation, were different within the same laminate. 

There were also, however, some common features among 
the stress distributions in different laminates. For ex- 
ample, the position of the maximum radial stress within any 
individual lamina was at the point where fibers are perpen- 
dicular to the hole edge, and* the maximum tangential stress 
was always located at the point where fibers are parallel to 
the hole edge. Stress concentration factors are also changed 
by varying the level of applied load. This is indicative of 
the non-linearity of the problem, since stress concentration 

if 

factors remain unchanged in linear cases. (See [1] for de- 
tails about these discussions.) 


if 

Numbers in brackets in this section refer to the references 
which are listed on page 83. 


81 


c. Progress During Report Period 

It would appear that the concept used for a single 
fastener is applicable to multi-fastener cases. An analysis 
model for a multi- fastener joint has been established as 
follows: The joint model is pictured as having a single row 

of equally spaced fastener holes infinite in number. This 
row of holes is parallel to the end of the joined plate 
(Figure IV-A-2-b) . As shown in this figure, the mathemati- 
cal model (see enlarged section) can be established by iso- 
lating a representative fastener hole. The plate width, W, 
of this model is taken as equal to the fastener spacing, s. 

The boundary conditions for the model of Figure IV-A-2-b 
are the same as for the single fastener hole case (Figure IV- 
A-2-a) , except for those along the boundaries where x = ±W/2. 
On these particular boundaries it is clear in Figure IV-A- 
2-b that displacements in the x~direction must be equal to 
zero, whereas the normal stresses are undetermined. On these 
boundaries in the single-fastener joint model of Figure IV- 
C-2-a just the reverse was true. Note that, in both cases, 
however, shear forces along these boundaries must be zero. 

In order to make comparisons between the stress behaviors 
of these two cases, we plan to set all the geometrical and 
laminate parameters of this model to be exactly the same as 
those of the previous single-fastener model. In addition, 
since there are no geometrical changes, we will use the iden- 
tical FEM mesh which was used for the single-hole case. The 
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applied loads will be chosen such that the load transmitted 

per pin will be 4,000 pounds and 6,L'00 pounds for B-type 

([±45®/0“] ) and 4,000 pounds and 5,760 pounds for P-type 

ns 

( [±30®/0®] laminates as before. 

XlS 

Rather than repeat all three analyses, however, we will 
use the minimum energy potential technique exclusively in 
the analysis of the multi-fastener joint. 

d. Plans for Upcoming Period 
We plan to complete the analyses which have been out- 
lined in this rep^art and to compare those results against 
those from the single-hole case. 

Wo will also plan to explore the application of the 
technique developed in this research to interference-fit 
cases . 


e. References 


1. Kim, W., "Stress Analysis Methods for Clearance-Pit 
Mechanical Joints in Laminated Composites", Ph.D. 
Dissertation, Rensselaer Polytechnic Institute, Troy, 
New york, 1982. 


3, Compact Lug Design 

Senior Investigator: D. B. Goetschel 

a. Introduction 

A critical aspect of many aerospace structural elements 
is the load transfer that takes place between the connecting 
lugs at the ends of the structure and the portions of the 
structure wherein the loads are well-distributed if not uni- 
form. An example is provided by the engine drag strut of 
■Che L-1011 aircraft (see Figure IV-A-3-a) . Such lugs are 
highly loaded and have a very complex stress state {see Fig- 
ure IV-A-3-b) . Further, since they must mate with connecting 
parts, these lugs are usually designed within rather strin- 
gent dimensional envelope constraints. As compared to lugs 
which (as in the case of the L-1011 strut) are made from 200 
ksi steel, meeting the geometric constraints, even with a 
composite structural design making maximum use of unidirec- 
tional graph! te -epoxy , could prove to be a difficult task. 
This research is intended to solve the load-volume problem 
in favor of composites rather than either reverting to de- 
signs using 200 ksi steel lugs attached to a composite strut 
or forcing the dimensional constraints to be relaxed with a 
redesign of the mating parts. 

b. Status 

Several initial test lugs using various ply layups were 
manufactured and tested. The results from these tests 


preceding page 


blank not FiEMEQ 




Dir^NSiONS 



MAXIMUM TENSILE LOAD 
MAXIMUM COMPRESSIVE LOAD 
MAXIMUM LATERAL LOAD - 


INCHES 


Figure IV-A— 3-b. Lug Experiment: Full-Scale L-1011 
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suggested some design modifications; namely, using bushings 
bonded inside the hole and slots at the sides of the hole. 

A second round of exploratory tests were performed to evalu- 
ate these design modifications. Some improvements could be 
seen, but the large scatter in the results made definite 
conclusions impossible. 

c. Progress During Report Period 

Graduate student Matt Cackett began work on the lugs in 
June. Initial testing was done using small quasi-isotropic, 
[0/45/-45/90] , strips loaded through a . 5-inch diameter pin 
hole. These initial specimens were not scaled to the Lock- 
heed dimensions; instead, available cured graphite plate 
material was used. Two slotted and one unslotted specimens 
were raade (see Figure IV-A-3-c) . 

An attempt was made to deter'" 'ne the stress cciicsntra- 
tion at the side of the hole experimentally, using Strass- 
Kote strain lacquer. However, it was not possible to accur- 
ately observe lacquer cracks. The test results are listed 
in Table IV *A-3-a. Note that the slotted specimen failed at 
a 22% higher load. It also failed in the bearing mode in- 
stead of net tension. 

Another attempt to determine the stress concentration 
at the edge of the hole was made using quarter-scale lugs 
with and without slots with a ply layup of [0/45/-45] (see 
Figure IV-A-3-d) . Photoelastic coatings were mounted on 
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TABLE IV-A“3-a 

STRENGTH COMPARISON OP VARIOUS LUG ARRANGEMENTS 


Specimen 

Number 

1 

2 

3 




unslotted 

unslotted 

slotted 


Failure Load 
(lbs.) 

4370 

4265 

5475 


these specimens. This technique was used effectively by 
Wonsub Kim for determining stress distributions in bolt- 
loaded plates. Residual stresses from tight fitting bolt 
connections and the small coating area in the region of in- 
terest contributed to inconclusive results from this method. 

It was decided to make a 12-specimen comparison between 

slotted and unslotted lugs. To simplify T:he interpretation 

of differences , both types should be made to fail in net 

tension. Thus, we desired a layup whose strength approached 

that of samples used during the previous reporting period, 

which did fail in net tension. A layup of [0/45/-45] was 

s 

chosen for this slotted versus unslouted comparison. Twelve 
specimens of 5208/1048 AIE graphite/epoxy were cut from the 
same plate. The cured plate thickness was found to vary 
from .454 inches near the center to .397 inches near the 
edges. This appeared to be due to using too thin (.125 in- 
ches thick) an aluminum plate on top of the laminate during 
autoclave cure. The aluminum deformed under the 100 psi 
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autoclave pressure resulting in uneven resin bleed-out. 

I 

Following removal of the plate from the autoclave, it was 
noticed that the bleeder material near the center of the 
plate was not saturated with resin. To ensure uniform 
clamping in the testing machine, one face of the lug was 
ground parallel to the other face over 1,5 inches in the 
gripped region only. 

Machining the specimens to size was done first using a 
carbide band saw to rough cut the lugs. Then the bolt hole 
was drill'ic' using successively larger high speed steel bits 
and, finally, a .5 -inch carbide tipped bit. Drilling was 
done at a feed rate of about .5 inches per hour. The slow 
feed rate was used to minimize hole surface flaws and break- 
out damage. No breakout was detected at this speed, and the 
hole surface appeared smooth. The outside surface was then 
ground to size using a heat treated steel template pinned to 
the lug as a guide. The slots were cut using a .03- inch 
thick diamond wheel mounted on a routing machine. Cuts were 
made from each side of the lug to maintain symmetry. 

The twelve specimens were tested on an Instron by clamp- 
ing one end and applying a load through the hole with a hard- 
ened pin. The test results were unexpected, in that the 
slotted lugs failed at an average load only three percent 
higher than the unslotted lugs. Finite element analysis re- 
sults had predicted that the tensile strength should be about 
twenty percent higher, based on peak tensile stresses. Three 
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points are worth noting, however. First, the results were 
reasonably consistent, as shown by the standard deviations 
listed in Table IV-A-3-b. The hxgher deviations for the 
slotted slugs may be attributed to the extra manufacturing 
process of slotting which was difficult to duplicate from 
piece to piece. Second, although the failure loads were 
quite similar, the failure modes were not quite the same. 
Both types failed in a combined tensile and shear tear-out. 


TABLE IV-A-3-b 

TEST RESULTS OP UNSLOTTED AND SLOTTED LUGS 


Unslotted Lugs 

Slotted 

Lugs 

Specimen 

Number 

Failure 
Load (lbs.) 

Specimen 

Number* 

Failure 
Load (lbs 

1 

12,125 

2 

13,250 

3 

12,220 

4 

11,210 

6 

12,970 

5 

12,715 

8 

12,715 

7 

12,360 

9 

11,480 

10 

14,120 

11 

12,830 

12 

12,125 

Average 

« 12,390 

Average 

= 12,630 

Standard 

Deviation 

= 559 

Standard 

Deviation 

= 99G 
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but the angles of the failure planes were different^ as 
shown in Figure IV-’A-3-e. This indicates that the unslotted 
specimens are tending more towards net tension failure than 
are the slotted ones. Thus, slots must reduce the tensile 
stresses by some palpable amount. Third, two lugs actually 
surpassed the Lockheed criteria (13,000 pounds~scaled) . One 
slotted lug even reach 14,000 pounds. The average failure 
loads for both types were within five percent of the design 
load. 

In order to isolate the effect of slots, twelve new 
specimens were made of a quasi-isotropic layup ([0/45/-45/ 

90] 2 q) with square extended heads (see Figure IV-A-3-f ) . 

This modification is a further attempt to reduce the possi- 
bility of shear failure by increasing the area of the shear 
plane. It was found, however, that these 16-ply specimens 
are susceptible to bearing failure at very low loads. They 
re currently being reduced in width to force net tension 
failure . 

The use of peak stress in failure criteria for the bolt- 
loaded lug appears to be a conservative approach. An alter- 
native failure criterion for bolt-loaded holes has been pro- 

[ 1 ] * 

posed by Agarwal . His criterion assumes failure will 
occur when the average stress over a certain distance is e- 
qual to the laminate strength . This approach also allows 

* ! 

Numbers in brackets in this section refer to the references 
which are listed on page 100. 
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Figure IV-A-3-e. Comparison of Angle of Failure Planes 

for Unslotted and Slotted Lugs 
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for tho three modes of failure; not tension# shear and bear- 
ing. This theory is now being applied to the aforementioned 
16-ply lugs using a finite element analysis. 

A capstrip (consisting of only 0°-plies wrapped around 
the outer edge of the lug) modification is now being inves- 
tigated. Instead of using a capstrip across the entire net 
tensile area, it was suggested by undergraduate student Tom 
DeMint to bond a partial capstrip to the lug (see Figure IV- 
A-3-g) . The advantage of this geometry is that the partial 
capstrip would increase the strength of the tensile cross 
section while retaining the bearing properties of the inner 
laminate. This design has the potential to be optimized by 
varying the partial capstrip width until the stress concen-- 
tration at the interface causes the capstrip to fail at the 
same instant the inner laminate fails in tension (unless, of 
course, some other failure mode has become active) . The in- 
creased stiffness of the capstrip should reduce the stress 
concentration at the edge of the hole, allowing the inner 
laminate to fail at a higher average load level. Initial 
testing is being done with a partial capstrip of .10 inch 
thickness. 

A lug with a partial capstrip was tested earlier in this 
project by Christoph Muser*- . A steel bushing was used to 
take the bearing load, and the capstrip extended from this 
bushing, as shown in Figure IV-A-3-gC. This lug failed in 
tension at a net tensile stress of 92,857 psi. 
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Following up on these promising results, a partial cap- 
strip of .10 inches was tested, failing at 12,750 pounds or 
36,890 psi. The capstrip sheared from the inner laminate in 
the bf'nd area along the sides of the lug away from the end. 
By improving the bond, the capstrip should produce good re- 
sults . 


d. Plans for Up c oming Period 

Plans for the coming period include improving the cap- 
strip design by strengthening the bond and optimizing the 
thickness. Also, a softening strip design will be investi- 
gated with use of Kevlar or glass inlays near the hole, in 
place of 0 ° -graphite-plies . Undergraduate student Bob Benson 
is pursuing a similar design using metal inlays. 

An accurate finite element model will be completed, 
either by using a cosine normal load distribution or by 
modeling the pin itself. Results from this will be used in 
conjunction with Agarwa]. 's failure criterion to determine 
the criterion's usefullneiss in predicting failure loads of 
compact lugs . 

'31 

A paper by T. A. Col Lings'- ■’ indicates that a layup of 
[0/45/0/-45] may be more effective in all modes of failure 
than the [0/45/-45] layup used in the slotted versus unslot- 
ted test. This n.ew layup also has been shown to have a 

. f41 

twelve percent higher peak stress concentration'- ■' than the 

[0/45/-45]. This would highlight the effect of slotting the 


lugs and may be pursued for another quarter scale test. It 
appears that the effectiveness of slotting may be more of a 
function of layup than was originally expected. The focus 
of continued work will be to quantify the effect of slotting 
on a reasonable layup and to evolve a design that will satis- 
fy a high load lag requirement (as typified by the Lockheed 
criteria) . 


e , References 

1. Agarwal, B, L., "Static Strength Prediction of Bolted 
Joints in Composite Materials", AIAA Journal, Vol. 18, 
No. 11, November 1980, pp. 1371-1375. 

2. "Composite Structural Materials", RPI Semi-Annual Pro- 
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April '79, pp. 36-41. 
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4. failure Criteria for Stress Concentrations 

Senior Invescigator: D. B. Goetschel 

a. Introduction 

It has been observed that when holes are drilled in ad- 
vanced composite plates whose dimensions are very large com- 
pared to the hole diameter, the failure load varies with 
hole sir.e, even though the elastic solution for peak stress 
next to the hole is independent of hole size. This shows 
that there is some mechanism other than peak stress contri- 
buting to failure when stress concentrations exist in struc- 
tures made from advanced composites. Some reduction of the 
peak stress due to non-linear material response would be ex- 
pected, but large variations have been found in advanced 
composites which are too brittle for this to be a reasonable 
explanation. Since lugs contain large stress concentrations, 
this behavior must be taken into account when they are anal- 
yzed. 

The work described here has been performed by Fola 
Ogunlari in research conducted in partial satisfaction for 
a Masters Thesis. 


b. Status 

The most attractive and general way of dealing with 
this phenomenon (i.e., where stress concentrations affect 
failure load) appears to be to use an average stress 
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r T_ 1 ^ 

criterion as doacribod by Whitney and Wuismer ^ ^ and by 
The essential feature of this criterion is use of 
an average stress calculated over a characteristic distance 
along a stress gradient. When the average stress, calcu- 
lated in this way, exceeds the failure stress for the mate- 
rial in a uniform stress field, it is postulated that the 
specimen will fail. The rational for this is as follov;s ; 

(a) brittle failure of a body under a given stress field is 
generally attributed to the existence of inherent flaws of 
various dimensions distributed throughout the body ; (b) 
Woibull statistics estimate the effect that those flaw 
distributions in a given material have on the failure 
strength of particular specimens and (c) it follows that the 
sharper the stress concentration, the less likely it is that 
a critical flaw will be encountered by the peak stresses, 
and thus, the higher the peak stress can go before failure 
occurs. The use of the average stress defined above is sim- 
ply an approximate method that appears to have some relia- 

. . [51 

bxlity in quantifying this behavior. Agarwal * applied 

this criterion to a particular bolted joint design and found 

that it accurately predicted static strength. 


Numbers in brackets in this section refer to the references 
wnich are listed on page 106. 
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c. Progress During Report Period 

Before the VJhitney/Nuisnier - Wu average stress criter- 
ion can be applied to the structural design process, it 
should first be tested for greater generality than has been 
evident in the literature to date. Accordingly, exact elas- 
tic solutions were sought and found for various types of 
(empty) holes and notches in isotropic elastic plates under 
load. From among these, solutions were selected for hole 
geometries that could be conviently manufactured and which 
provide a full range of stress concentration intensity. 
Figure IV-A-4-a shows the configurations selected. As the 
diameters of the round holes increase relative to the width 
of the test strips, the stress fields approach uniformity. 

In the other extreme, the square ended notches provide a 
stress singularity. 

The specimens were manufactured from 1048 AIE graphite/ 

epoxy using a layup of [0/90/45/-451 . To determine the un- 

s 

notched strength and the characteristic dimension for inte- 
gration, a series of tests were conducted on unnotched and 
round-holed specimens. A minimum root-mean-aquare error 
analysis resulted in a characteristic dimension of .195 in- 
ches. Table IV-A-4-a shows the results of the final test 
sequence. As can be seen, the average stress criterion was 
found to be ^’’ery reliable over a wide range of stress con- 
centration levels. 
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TABLE IV-A-4-a 

SUMMARY OP AVERAGE STRESS AND PEAK RESULTS 


Hole and 
Crack 
Geometry 

Test 

Failure 

Load 

(pounds) 

Average 

Stress 

Prediction 

(load/ 

% error) 

Peak 

Stress 

Prediction 

(load/ 

% error) 

6/32 dia. 

3450 

3712 / 

7.6 

190S / -44.7 

9/32 dia. 

3038 

3183 / 

4.8 

1781 / -41.4 

15/32 dia. 

2042 

2393 / 

17.2 

1392 / -31.8 

Double edge 

2760 

2884 / 

4.5 

0 /-lOO.O 

Center crack 

2633 

2619 / 

-.5 

0 /-lOO.O 


Plans for Upcoming Period 

The dependence on plate thickness of* the characteristic 
integration length will be the next subject for investiga- 
tion. If the rational for using an average stress criterion 
is based on the probability of a critical sized flaw occur- 
ring in the area of the largest stresses, then thicker plates 
should increase that probability and shorter characteristic 
lengths should be used with thicker plates, in recognition 
of the fact that failures will occur sooner relative to an 
uniform stress material allowable. 

If this is successful, the method will be exercized in 
analyzing lugs made using advanced composites. 
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IV-B SCALING EFFECTS IN TESTING COMPOSITE STRUCTURES 

Senior Investigator: D. B. Goetschel 

1, Introduction 

It can bo extremely valuable in terms of both time and 
monay to bo able to tost small scale models when designing 
largo structures. This is especially true for structures 
fabricated from composite materials such as graphite/opoxy , 
since the material itself is relatively expensive. On the 
other hand, the knowledge required to interpret results from 
small scale models, accurately and with confidence, is essen- 
tially nonexistent, mostly because the effects of the unique 
characteristics of composite materials are only partially 
understood. 


2. Status 

Undergraduate student Scott Baxter conducted a prelimi- 
nary set of tests for a series of samples with different 
scales, loaded in simple static tension. Good correlation 
was found for strength variation with scale between experi- 
ment and the theory that was developed for the limited num- 
ber of tests that were run. The predictions were based on 
free edge and Weibull statistical effects . 

3. Progress During Report Period 
Wherever statistics are involved, there will be a ques- 
tion as to the adequacy of the number of samples tested. 
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The preliminary testa, reported earlier, wore considered too 
few. Work has been carried out to conduct a final set of 
verification tests for simple static tension, which will 
broaden our statistical base. 

4. Plans for Upcomming Period 
A final set of tests on simple static tension will be 
completed and compared witl the developed theory. Work will 
be started on the scaling of bearing failures at pin-loaded 
holes. Exploratory tests are also being planned to investi- 
gate in a broader sense the scaling of experiments on joints 
in composite materials. 
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IV«-C END EFFECT DECAY IN PRISMATIC MEMBERS 

Senior Investigator: D, B, Gootschol 

1 . Introduction 

Saint Venant’s Principle states that all statically 
equivalent systems of loads on a body will produce nearly 
identical stress fields in regions that are remote from the 
loaded area. The usual engineering approximation of neglec- 
ting Saint Venant's end effects at distances of about one 
width from the ends is justified for isotropic members of 
compact cross section. For composites like graphite/epoxy, 
however, or for open thin-walled cross sections, these ap- 
proximations can be substantially in error. It is often 
important to be able to analytically quantify these effects, 
in general, for any material system and cross sectional geo- 
metry. It is particularly useful to be able to provide this 
type of analysis for long, uniform struts such as are found 
in structures like the L-1011 aircraft, wing-mounted engine 
drag strut. Much shorter specimens could be used for test- 
ing, if such analyses could define minimum lengths above 
which stress distributions are not altered. Work by graduate 
student Tsay-hsin Hu to accomplish this analysis is des- 
cribed here, but first some background information will be 
given . 
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2 . Status 

n 1 * r 21 

Toupin*''^'' and Knowles^ •* independently prosentuse. theo- 
rems for upper bound estimates of the strain energy 

V(x) => V(0) *exp(-2Xx) 

where: X = characteristic decay rate of end effects, 

V(0) = total strain energy and 

V(x) 5= strain energy in body beyond x. 

Then, because of the quadratic nature of strain energy in 
terms of the mechanical variables, the displacements, strains 
and stresses relationships are, respectively, 

u^(x) <_ Uj^ (0) *exp(-Xx) 

Eij(x) < (0) *exp(-Xx) 

^ij (x) 5 (0) *exp(-Xx) 

Choi and Horgan^^^^'^'^^ evaluated X* for some specific 
cases by employing the Airy stress function c|) 

4)(x,y) = F(y) *exp(-Xx) 

This solution form enables the mathematics to be reduced to 
an eigenvalue problem involving a fourth order ordinary dif- 
ferential equation in F (y) . Since all the exponents are 
negative, the eigenvalue with the smallest real part charac- 
terizes the dominant exponential decay rate. 

f 71 

Dong and Goetschel developed a method for finding 


Numbers in brackets in this section refer to the references 
which are listed on page 117. 


Ill 


the decay rate of edge effects in a plate composed of an ar- 
bitrary number of anisotropic layers. Instead of the Airy 
stress function, the in-plane displacements were adopted as 
the primary independent variables. One-dimensional finite 
element polynomial interpolation functions were used to model 
the behavior through the thickness of the plate. The com- 
puter program "Plate" was the result of this approach to 
quantifying edge effects in serai-infinite plates. "Plate" 
has a serious limitation, however, in that It can only deal 
with plates with orthotropic layers. This is because it 
lacks displacement degrees of freedom in the direction paral- 
lel to the free edge of the plate. 

A program called "Prism" was partially developed subse- 
quent to "Plate's" being established? its purpose was to 
quantify end effects in prismatic members, as shown in Fig- 
ure IV-C-1. By restricting the degrees of freedom it can 
handle, "Prism" is also to model general anisotropic layered 
plates. Development of this program reached the point of 
assemblying the governing system matrices , but it was not 
completed to the point of solving the resultant eigenvalue 
problem. 

"Prism" uses a semi-analytic finite element method to 
quantify the end effects. Two-dimensional finite element 
polynomial interpolation is used to model displacement vari- 
ations across the member cross section. Then an exponential 
decay of displacement is assumed along the length of the 
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mombor. The rosulting oquationo form a aocond order eigen- 
problem, whose roots ere the characteristic decay rates and 
whoso corresponding eigenvectors are the diaplacomont distri- 
butions across the cross section, The complete displacement 
field is then obtained by multiplying the cross sectional 
variation by the axial exponential decay. The stress field 
is easily obtained from the displacement field as a final 
stop. 

The semi-analytic finite element formulation assumes a 
diaplacomont field of the form 

u(x,y,z) “ p(y,z) *oxp(-Xx) 

whore: p(y» 2 ) » polynomial interpolation across the 

cross section and 

exp (-Ax) a exponontial decay interpolation along 
the length. 

A typical four-nodod finite element ’is shown in Figure 
IV-C-2. 

The usual finite element forco-diaplacoment matrix 
equation is 

tk] {u} a (f) 

whore: [k] >=' stiffness matrix, 

Cu) a vector of displacement degrees of freedom and 
(f) ° vector of force degrees of freedom. 

The stiffness matrix is composed of three parts; 


[k] [kp] + [kg] [k^l 




the small arrows marked from I to \Z are 
the finite element degrees of freedom 


Figure IV-C-2. 

Configuration of a Typical Four-Noded Finite Element 
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Here [k ] is calculated by integrating tractions over the 
end of the element, [k„] is calculated by integrating trac- 

O 

tions over the sides (including top and bottom) of the ele- 
ment and [k^] is calculated by integrating body forces over 
the volume of the element. The activation of any one of the 
displacement degrees of freedom defines a displacement field 
throughout the element. This displacement field can he used 
to find the corresponding stress field. The stress field, 
in turn, defines the tractions that are integrated over the 
surface of the element. The derivatives of the stress field, 
when inserted into the differential equilibrium equation, de- 
fine unbalanced body forces that are integrated over the 
volume of the element. The displacement interpolation func- 
tions are used as weighting functions when performing the 
various integrations. 

The problem under consideration requires that there not 
be any tractions on the sides or body forces within the in- 
terior of the element, therefore, 

([kg] + [k^l){u} = {0} 

The differentiations and integrations cause [kg] and [k^] to 
be made up of three different matrices with different powers 
of X 

[kg] + [k^l = XEk^] + [k^3 + ^[k^l 
Thus , the eigenproblem to be solved can be written as 
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(XMk2l + X[k^] + [kQ]){u} = {0} 

The solution for eigenvalues is complicated by the fact that, 
although [k^] and [k^] are symmetric, [k^^] is anti-symmetric. 

3. Progress During Report Period 

As a first step, "Plate" was transferred successfully 
into the RPI computer system. This was to provide a basis 
of comparison for the newly developed program "Prism" . 

"Prism" can be used to model plate problems, so intermediate 
results from the two programs during various steps of the 
solution process can be compared as an aid in debugging 
"Prism". "Prism" is now ge‘, aerating the same system matrices 
as "Plate". 

Despite the fact that the system matrices are the same, 
the new program was not getting answers in agreement with the 
old one. It was discovered that some of the system matrices 
are singular and the frequency shift technique is needed to 
overcome this problem. "Plate" uses a solution algorithm 
that performs frequency shifts automatically, however, this 
solution algorithm is very unreliable and frequently produces 
spurious results. It is not being used in "Prism". 

Although the singularity problem would seem to have been 
solved, consistent results are still not being produced by 
"Prism" . It appears that the problem is ill-conditioned for 
the solution methods that have been tried to date, because 
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some correct results have been produced, and the accuracy 
of the results decays as the number of degrees of freedom Is 
Increased. 


4. Plans for Upcoming Period 
VJork will be continued on the solution algorithm for 
"Prism" until consistent results can be obtained. Once this 
is accomplished, the program will be verified for prismatic 
members by comparing it to available results from exact 
analytic solutions . The program will 
then be used to investigate several problems of interest. 

In layered composite plates, the variation of edge effects 
as ply orientation and stacking sequence are varied will be 
examined. For prismatic members, the end effects of cross- 
sectional geometry (i.e., I-beams versus round tubes versus 
rectangular beams, etc.) will be investigated. 
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V"A INITIAL SAILPLANE PROJECT: THE RP-»1 

Senior Investigators; F. P. Bundy 

R. J. Diefendorf 
H . Hager up 
H, Scarton 


1. Status 

Recognizing that the resins in composite materials can 
change properties somewhat with aging and exposure to tem- 
perature/moisture cycles, the wing/fuselage ensemble struc- 
ture of the RP-1 glider has been static tested about once a 
year since its initial fabrication in '79- '80. The last 
such test was done at RPI on August 5, 1981, preparatory to 
winch launch flight tests at Schenectady County Airport. 

That test consisted of simple bending using a distributed 
load up to a maximum of 590 pounds per wing, which corre- 
sponds to a little over 4 g's for average pilot weight. 

2. Progress During Report Period 
After the failure of the wing beam "carry-through” 
structure of the RP-2 glider in November 1981, the corre- 
sponding parts of the RP-1 glider were studied with an eye 
to making it stronger and more fail-proof. Two changes were 
made in the carry-through structure: (i) the female open 

box member was wrapped with an additional layer of 45“/45“ 
Kevlar fiber and bonded with resin and (ii) the aluminum pin 
which goes through both wing stubs at the midpoint of the 
carry-through region was replaced with a stronger steel pin. 
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and the lands on the beam webs on which the pin bears were 
made thicker. The latter change was made to insuto that the 
pin assembly could absorb any tension, generated from tenden- 
cies to pull the wing beams apart, such as occur with wing 
fore-aft bending. 

in the present tests, the wing-fuselage assembly wati 
subjected to simple bending, as in previous years but, in 
addition, to a bending-torsion test, which would simulate 
the stresses in a high speed pull-out. It was in this mode 
that the Rp-2 wing beam carry-through failed (see Section 
V-B) . 

As a safety measure to prevent catastrophic failure of 
the whole structure in case of a failure of one part during 
the test, the wings each had two "safety stops' under them - 
one at about half span and the other a little inboard of the 
wing tip. As the wing loading and corresponding deflections 
increased, the "stop~supports" were adjusted so that the wing 
and load could fall no more than i|" to i" before being ade- 
quately stopped and supported. 

In addition to monitoring bending and tv/isting deflec- 
tions, the wing root regions were monitored by attached mi- 
crophones for any "acoustic events" associated with the stress- 
ing of the structure. The microphones were in contact with 
the wing skin on the top and bottom sides of each wing near 
the root ribs and as close as possible to the carbon fiber 
capstrips of the wing beams. This AE monitoring was done by 
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Professor ll. Scar ton and graduate student G, Dobal, using a 
four-channel physical acoustics model 3400 AE analyzer sys- 
tem. 

The other changes on the glider made since the previous 
strength tests (i.e., closing in the bottom of the cockpit 
area with foam ribs and belly skin, and streamlining the 
skid and wheel) had no effect on the main structural ele- 
ments being tost and, although they were in place at the 
time, were irrelevant to those tests. 

The RP-1 was assembled for testing and placed in an in- 
verted position, with the fuselage booms resting on a sup- 
port about 36 inches above the floor. Deflection measuring 
stations were established at inboard, mid and outboard points 
on the leading and trailing edges of each wing, making a 
total of twelve stations (as shown in Figure V-A-1) . At 
each station a vertical steel rod, supported by a cast iron 
base on the floor, held a vertical steel metric scale at 
roughly the level of the wing's leading and trailing edges. 
Crossed with the vertical scale was a horizontal scale at- 
tached to the wing (i.e., at that same spanwise station). 

This arrangement made it possible to read both vertical and 
horizontal displacements as the wing was loaded. 

In the simple bending tests, the fuselage was held in a 
horizontal position and the ten-pound sand bag weights were 
laid on, one per foot of span, directly over the wing beam. 

In the bending/torsion tests, the fuselage was positioned on 
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NOTE; Sand bags were loaded synunefcrically and simultaneously 
left and right, starting at inboard end. 
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Figure V-A-1. 

Static Loading Tests of RP-1 Glider Wing for 
Bending, Torsion and Fore-Aft Deflections 
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a ten degree slant (nose low, tail high, inverted) and the 
sandbag weights placed on a line 14 cm aft of the wing beam. 
Thus, both torsion and fore-aft bending were induced, in ad- 
dition to vertical bending. These two static tests were 
survived very successfully. The deflections at the differ- 
ent simple bending loadings are plotted in Figure V-A-2. 

These deflections compare closely with those obtained in the 
test of August 5, 1981. Upon unloading, the deflection 
readings came back to within 0,1 to 0,3 cm of the initial 
values, indicating that the structure had suffered no per- 
manent deformation. 

The results of the bending/torsion test at the 4 g 
loading are presented graphically in Figure IV-V-3. Differ- 
ences in trailing edge deflections relative to the leading 
edge divided by the chord distance gave the angles of twist 
shown in this figure. By usual standards, this is a very 
stiff wing. 

A major concern, because of the RP-2 failure, was how 
much the wings deflect forward during the bending/torsion 
test mode. The horizontal scale readings of deflections at 
the six stations at the leading edge had to be combined with 
vertical scale readings, since the bending axis of the wing 
was tilted ten degrees to the horizontal. Thus, if there 
were no deformation of the wing in its own plane, the trans- 
verse deflections would cause the leading edge to move aft 
as measured by the horizontal scales. This proper combination 
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of these measurements led to an estimate of total forward 
deformation of the wing at the outboard end of 36 mm. To 
ascertain the degree to which carry-through structure might 
be carrying fore-aft bending moments, the aft wing pin was 
watched during the loading and unloading. If there was any 
pull-out of the pin from its socket, it appeared to be less 
than a fraction of a millimeter. 

A substantial amount of acoustic emission data was 
taken during the test. The most salient of the output oc- 
curred during the bending/torsion portion. The four sensors 
located near the points of maximum stress (wing spar center 
section carry- through) recorded AE events of significant 
amplitude only at loadings between three and four g's. The 
counts versus time plot presented in Figure V-A-4 showed 
some indications of being critically active, according to 
the ASTM E569 criterion. Figure V-A-5 reveals seven large 
amplitude events above 76 dB; none occurred before the three 
g load point. A Hsu-Neilson lead-break calibration directly 
opposite the sensor on the root rib (0.5 mm lead) produced 
an amplitude of 86 dB. This same lead break calibration was 
not resolvable at a distance of two feet from the sensor. We 
thus conclude that all data in Figures V-A-4 and -5 is pro- 
duced by local activity. In the absence of definitive indi- 
cations, the elastic behavior of the structure had led us to 
believe that these relatively strong AE events were more 
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likely due to relative stick-slip motion of mating (but not 
bonded) surfaces rather than breakage of stressed fibers. 

Following the static tests, winch-launched flight tests 
were carried out on July 23 and 30, 1982 at the Saratoga 
County Airport for measurements of sink rate versus flight 
speed at a gross weight of 340 pounds (F, Bundy, pilot) . In 
some of these flights, the air was stable enough to yield 
sink rate data of reasonable reliability over the speed 
range from 25 to 45 mph. The effect of aerodynamic fairings 
was to lower the minimum sink rate to about 200 FPM (from 
250 FPM in the unfaired configuration) and to increase the 
optimum glide slope from about 10 or 11:1 to 14:1. Flight 
tests were also made at a gross weight of 380 pounds (M. 
McCarron, pilot) on September 17 and 24, 1982, but, except 
for a 30 mph run and a 40 mph run, the air was convecting 
too much to yield meaningful sink rate data. Two rather 
creditable soaring flights were accomplished; the first, on 
July 23rd, was a 2600-foot, 20 minute thermal soaring flight, 
and the second, on September 24th, reached 1500 feet, with 
an endurance of 19 minutes under overcast conditions. 

Tests were conducted on five identical elements of two 
different kinds , typical of structural elements used on the 
RP-1 and hidden from view, in the Pall of *81 and five more 
in September '82. The two joint types are shown in Figure 
V-A-6. These tests are intended to provide the first indi- 
cation of possible strength degradation that, with time. 
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Type A (open foam ) 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Gloss (2 sides) 

V-40/828 with 
Cotton Flux 

1/8“ Red Foam 
(open face) (hi density) 

Kevlar ( 1 side) 


Type B (closed foam) 



Gloss ( 2 sides) 


Kevlar ( 2 sides) 


V-40/828 with 
Cotton Flux 


Figure V-A-6 . Two Types of Structural Joints Used in the 
RP-1 (lightly-loaded, low-cost composite 
structures) 
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might be encountered with the RP-1, All elements are stored 
under the same conditions as the aircraft. Comparative re- 
sults of the tests in '81 and '82 are shown in Figures V-A-7 
through 9. It appears that degradation with time is well 
within the spread in strength due to fabrication or mate- 
rials quality. No further conclusions can be drawn with 
confidence at this time. 

One unexpected result from these tests is that the 
joint to a single sandwich skin is stronger than that to a 
double sandv/ich skin. lu appears that in the former case, 
the resin penetrates the foam and bonds to the single outer 
skin. This results in tensile strength greater than the 
"peel strength" between inner skin and foam core which is 
called into play in the latter case. 

3. Plans for Upcoming Period 

The RP-1 has now logged 51 flights and about four hours 
of flight time. 

In its present form, the machine is a reasonable glider, 
but it can be improved further by streamlining the upper 
pilot area (smooth bubble canopy) and by making the ailerons 
more control effective. Steps to accomplish these changes 
will be considered during the next period. 

Testing of typical structural elements will continue in 
the year ahead to see if any supnortable trends can be es- 
tablished as to how such elements degrade with time. 
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Figure V-A-7. Results of Tensile Tests to Failure of the Structural Joints Shown 

in Figure V-A-6 








ORIGINAL PAGE IS 
OF POOR QUALITY 




— Type A 

(averaged over time) 




Type B 

(averaged over time) 


O Type A 
O Type B 






V-B SECOND SAILPLANE PROJECT; THE RP~2 


Senior Investigators; 


F. Bundy 

R# J. Diefendorf 
H« Hagerup 
C^ Muser 


1 . Status 

As reported earlier/ the RP-2 sailplane was statically 
tested in October 1981. It successfully survived the 5.9 g 
bending test, but during the combined bending-torsion tests, 
the wing carry-through structure failed. This failure 
caused secondary damage to the outer wings, since they were 
only supported at the fuselage, and the wing tips fell to 
the floor. The damage was analyzed and several stress and 
deflection calculations performed, as a result of which it 
was coi'^irluded that the failure initiated as a result of the 
carry-through structure being loaded in fore-aft bending and 
torsion, which it was not intended to carry. Repair tech- 
niques were developed for the outer wings. The splices to 
be used to restore the structural integrity of the broken 
capstrips of the wing I-beam were successfully tested on a 
separate test beam. The wing carry-through structure was 
redesigned to provide for a better vertical shear transfer 
mechanism and to prevent inadvertenjj^y^ becoming a load-path 
for torsion and fore-aft bending, per se. 
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2, Progress During Report Period 

In the current reporting period, rocons -ruction pro- 
ceeded to the following point. The broken caps trips (at 
about one-half span) of both wings have now been repaired, 
by employing the splice technique developed during the pre- 
vious reporting period. One wing has been rebuilt, except 
fcr the carry-through portion of its spar, and a portion of 
its top wing skin remains to be installed. An existing test 
beam v;as taken apart in the carry-through region and flaws 
were introduced to simulate the damage in that area as it 
existed on the actual wing. The testbeam carry-through 
structure was tested up to the 5,9 g limit of the test fix- 
tures. The load was kept at its highest level for three 
minutes, and there was no evidence of damage of any kind. 

The new carry-through connection has been designed to with- 
stand three times that load. Following this test the deci- 
sion was made to install the new design of the actual wing 
of the RP-2. This wing connection was partially rebuilt at 
the end of the reporting period. 

The fuselage was also modified to accommodate the larger 
parts in the new wing-connection. Testing of the compres- 
sion strut that is to be mounted in the fuselage between the 
two shear-torsion pins, -which carry the load from wing root 
rib into the fuselage, has begun. This strut is the member 
which stiffens deflections associated with fore-aft wing 
bending and, thereby, helps to ’prevent the' wing-carry through 
from carrying such loads. 
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3. Plans for Upcoming Period 
The reconstruction of the wings will continue until 
the end of the year. This includes rebuilding the wing con~ 
nection to ailerons, mounting the flaps, replacing the ac- 
tuator mechanisms and splicing the wingskins. The steering 
mechanism in the cockpit also needs to be partially rebuilt. 
The compression strut will be finalized, follov/ing testing, 
and will be fabricated and installed. Finally, a test fix- 
ture will be designed and built to allow static testing of 
the completed airframe. 
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PART VI 

COMPUTER SOFTWARE DEVELOPMENTS 
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VI COMPUTER SOFTWARE DEVELOPMENTS 

Senior Investigator: M. S. Shephard 

1 . Introduction 

The objective of this portion of the project is to pro- 
vide advanced and specialized computer "tools" for the anal- 
ysis and design of composite structural elements and the 
study of structures-related phenomena in composite materials . 
The major thrust to date has been and will continue to be in 
the finite element area, with effort directed at implementing 
finite element analysis capabilities and developing inter- 
active graphic preprocessing and postprocessing capabilities. 

2 . Status 

The results achieved in this area of research during 
recent reporting periods have been used for the numerical 
analysis of moisture effects and the numerical investigation 
of the micromechanics of composite fracture discussed in 
Parts III-D and III-E of this report. In addition, the 
POPES* system, which is available to the composites program, 
continues to be developed through student class projects and 
projects supported by the Graphics Center's Industrial Asso- 
ciates Group. 


'ft 

POPES (People Oriented Finite Element Software) is RPI's 
finite element software sy^stem. 
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3. Proqrosa During Report Period 

Projects of interest to the composites program supported 
out of the Center for Interactive Computer Graphics and the 
Department of Civil Engineering are discussed in the follow- 
ing two subsections. 

g. Throe-Dimonaional Preprocessor, Based on Extrusion Capa - 

The extrusion preprocessor allows for the definition and 
meshing of three-dimensional objects by first defining com- 
ponents with planar faces and then extruding them in the 
third dimension and patching them to similarly-defined com- 
pononts ^ ^ . During the last reporting period, the capabil- 
ities of this program have been expanded to include more ox- 
tensive patching modes. In addition, the mesh generating 
capabilities are being improved. 

In addition to the simple face to face component patch- 
ing available previously, it is now possible to patch new 
components to planar sides. Also, the ability to do linear 
extrusion in which a face is extruded in a linear manner to 
a new face that is a scaled version of the original has been 
added. Figure VI-1 displays an example of the type of geo- 
metric object that can be generated witli these capabilities. 

The mesh generating capabilities have been extended to 
allow for the meshing of component faces in terms of any 


Numbers in brackets in this section refer to the references 
which are listed on page 151, 
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Figure VI-1 Example of Patched Components 
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number of subregions (Figure VI~2) which can then be extruded 
to form a solid mesh within that region (Figure VI-3) , The 
meshes for various compatible subregions can then be patched 
together. Figure VI-4 shows a two-component object fully 
meshed. 

b. Composite Plate and Shell Element 

As discussed in the previous progress report, an ele- 
ment was being established in our in-house finite element 
analysis program for the analysis of composite plates and 
shells. The element selected is a flat triangular element 
which results from the superposition of the DKT plate bend- 

r21 

ing element and a constant strain triangle for membrane 
behavior. 

The DKT plate bending element was first incorporated 

and tested. The results obtained for this case were very 

good. The next steps were to add the constant strain tri- 

(31 

angle and the laminated plate capabilities ■' . Both of 

these capabilities have been tested and good results ob- 
(41 

tained . The final steps in the element derivation include 

providing the bility to place this element in a general 

three-dimensional domain and add a corrective stiffness term 

to account for the rotational degree of freedom normal to 
(51 

the element 

To facilitate the testing of the shell element, a link 
is currently being created between the shell preprocessor 
and the analysis program. This link also allows for the 
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Figure VI-3 Three-Dimensional Mesh in Extruded Component 
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interactive alpha-numeric specification of problem loads, 
material propeirties and boundary conditions. The analysis 
results are also being linked to the three-dimensional dis- 
placement postprocessor. Initial shell analysis runs in- 
dicate satisfactory results, although they are somewhat less 
accurate than the bending results. This was expected, since 
the constant strain triangular representation of membrane 
behavior is only an approximation. It is very difficult, 
however, to develop combined elements for shells with better 
than linear membrane models, 

4, Plans for Upcoming Period 

During the next reporting period additional features 
will be added to the extrusion program and shell preprocessor 
{discussed in previous progress reports) . In addition, a 
final set of test problems will be run on the composite flat 
shell element added to the in-house finite element program. 

An investigation is also being initiated to investigate 
the use of modified-quadtree mesh generation capability 
developed at RPI in a finite element analysis program to 
automatically track down propagation. If such a procedure 
is successful, it could have a large impact on the efforts 
in the numerical investigation of micromechanical composite 
fracture discussed in Part III-E. 
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1 

1 

I 

> 

INTRODUCTION 

Technical meetings, on~ and off-campus, provide impor- 
tant opportunities for interchange of technical information. 

Because of the large number of composites meetings, a cen- 
tral catalog with all upcoming meetings is being maintained 

5 

and distributed periodically. In this way we help to assure 
that a Rensselaer faculty/staff member can participate in im- 
portant meetings. The calendar for this reporting period is 
shown in Table VII-1. Meetings attended by RPI composites 
program faculty/staff/students during the reporting period 
are shown in Table VII-2. Some meetings particularly rele- 
vant to composites , held on-campus with special speakers , 
are listed in Table VII-3. A list of composite-related vis- 
its to relevant organizations by RPI faculty/staff /students, 
with the purpose of each visit outlined, is presented in 
Table VI 1-4. 

• Two continuing education special courses , outgrowths of 
the composite materials and structures program, were pre- 
sented for graduate engineers in industry and government, 
for the third time: "Advanced Composite Materials and 

Structures" during the week of July 12-16 and "The Finite 
Element Method" during August 2-6, Each course lasted a week, 
and the level of the material was again planned to be particu- 
larly useful to managers of engineering structures activities 
who are involved in technical work but who may not have taken 


• *•» 
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TABLE VII -I 

CALENDAR OF COMPOSITES-RELATED MEETINGS 
(April 30, 1982 through September 30, 1982) 


1982 

5/4-7 38th Annual Forum and Technical Di3play, Anaheim, 

CA. “Sponsored by AHS." 

5/10-12 Structures, Structural Dynamics and Materials 

Conference, New Orleans, LA. "Sponsored by AIAA/ 
AHS/SDM/ASME . “ 

5/13-14 Workshop on Modeling, Analysis and Optimization 
Issues for Large Space Structures, Williamsburg, 

VA. "Sponsored by Air Porce/NASA." 

5/17-19 2nd International Very Large Vehicle Conference, 
Washington, DC. "Sponsored by AI7VA. " 

5/17-21 Composite Metals Workshop, Berkeley, CA. "Spon- 
sored by the University of California." 

5/19 Polymer Symposium, Urbana, IL. "Sponsored by the 

University of Illinois." 

5/20-21 13th Akron Polymer Conference on Fracture of Poly- 
mers and Composites Program, AJcron, OH. "Sponsored 
by ACS . " 

5/25-27 Annual Meeting and Technical Display , Baltimore, MD. 
"Sponsored by AIAA. " 

6/14-18 Gordon Conference on Chemical Characterization of 
Structural Polymers, Plymouth, NH 

6/16-rl7 Annual NCGA Conference, Anaheim, CA. 

6/21-25 9th U. S. Congress of Applied Mechanics, Ithaca, NY. 
"Sponsored by AIAA. " 

7/12-16 lUPAC, Amherst, MA. 

7/13-14 Symposium on Jointing in Fibre Reinforced Plastics, 
London, UK. "Sponsored by the Aeronautical Depart- 
ment, Imperial College and the Royal Aircraft Es- 
tablishment . " 

7/15-19 2nd Annual Computer Engineering Conference, San 
Diego, CA, "Sponsored by ASME." 

7/16-19 Symposium on Mechanics of Composite Materials, 
Blacksburg, VA. "Sponsored by lUTAM." 

7/22-27 13th Congress, Seattle, WA. "Sponsored by ICAS/AIAA. " 
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8/23-27 

8/30-9/2 

9/16-17 

9/20-24 

9/21-23 

9/24-30 


Gordon Conforenco on Interfaces, New Hampton, NH. 

lUPAC Polymer Processing Symposium, Athens, Greece, 

CFRC, Bonbanne, Prance. "Sponsored by the French 
Carbon Group." 

Carbon '82, London, UK, "Sponsored by the Society 
of the Chemical Industry." 

Symposium on Solid Mechanics, S. Yarmouth, MA. 
"Sponsored by the U. S, Army." 

Aerospace II Symposium, Williamsburg, VA. "Spon- 
sored by NAA, " 


9/27-10/2 33rd International Astro. Congress, Paris, Prance. 

"Sponsored by the French Aero, and Astro. Association." 
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TABLE VI 1-2 

COMPOSITES-RELATED TECHNICAL MEETINGS ATTENDED OFF-CAMPUS 
(April 30, 1982 through Sep f;.:jinber 30, 1982) 


1982 

5/13-14 Air Porce/NASA Workshop on Modeling, Analysis and 
Optimization Issues for Large Space Structures 
(Prof. Loewy) , Williamsburg, VA. 

5/19 University of Illinois Polymer Symposium (Prof. 
Sternstein) , Urbana, IL. 

Professor Sternstein gave invited lecture; 
"Deformation and Failure in High Perform- 
ance Composites". 

5/20-21 13th Akron Polymer Summit Conference on Fracture of 
Polymers and Composites Program (Prof. Sternstein) , 
Akron , OH . 

Professor Sternstein gave invited lecture; 
"Fracture in High Performance Laminates". 

5/25-27 AIAA Annual Meeting and Technical Display (Prof. 
Loewy) , Baltimore, MD. 

Professor Loewy chaired a session on cruise 
missies . 


6/14-18 Gordon Conference on Chemical Characterization of 
Structural Polymers (Prof. Sternstein) , Plymouth, 

NH. 

Professor Sternstein gave invited lecture; 

"The Influence of Resin Chemistry on Mechan- 
ical Properties" on June, 15 and was the 
morning discussion leader on June, 16. 

6/16-17 Annual NCGA Conference (Prof. Shephard) , Anaheim, CA. 

Professor Shephard presented; "Interactive 
Finite Element Programs" and "Finite Element 
Models from Geomebric Models". 

6/21-25 9th U. S. National Congress of Applied Mechanics 
(Prof. Brunelle and Dr. Muser) , Ithaca, NY. 

Professor Brunelle presented paper; "Eigen- 
value Similarity Roles for a Class of Rec- 
tangular Specially Orthotropic Laminated 
Plates" . 
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7/12-16 lUPAC (Profs. Diefendorf emd Sternstein ) , Amherest, 
IIA. 

Professor Diefendorf co-chaired a session 
on composites. 

Professor Sternstein gave lecture; "Visco 
Elastic Characterization of Neat Resins 
and Composites". 

8/16-19 lUTAM Symposium on Mechanics of Composite Materials 
(Dr. Muser ) , Blacksburg, VA. 

8/17-18 2nd Inter. .ational Computer Engineering Conference 
(Prof. Shephard) , San Diego, CA. 

Professor Shephard presented paper; "An 
Approach to Automatic Finite Element Mesh 
Generation" and, in the Product Design 
with CAD Systems Panel, presented; "Pre- 
paring Mechanical Engineering Students for 
the CAD Environment" . 


8/23-27 Gordon Conference on Interfaces (Prof. Diefendorf) , 
Nev/ Hampton , NH . 

Professor Diefendorf chaired a session on 
composite interfaces. 

8/30-9/2 lUPAC Polymer Processing Symposium (Prof. Stern- 
stein) , Athens, Greece. 

Professor Sternstein gave lecture; "Defor- 
mation and Failure of Composites". 

9/16«rl7 CFRC of the French Carbon Group (Prof. Diefendorf) , 
Bonbanne , France . 

Professor Diefendorf gave invited lecture; 
"Discotic Liquid Crystals". 

9/20-24 Carbon '82, Conference of the Society of the Chem- 
ical Industry (Prof. Diefendorf), London, UK. 


NAA Symposium, Aerospace II (Prof. Loewy) , Wil- 
liamsburg, VA. 

Professor Loewy chaired a session on aero- 
nautical research and development. 


9/24-30 
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TABLE VI 1-3 

COMPOSITES-RELATED MEETINGS / T ALKS HELD AT RPI 
(April 30, 1982 through September 30, 1982) 


Tooic 


A Survey of Advanced Com- 5/4 
posite Applications to 
Military Aircraft 

A University Program of 5/17 
Composite Sailplane De- 
velopment 

An Incremental Life Predic- 6/16 
tion Law for Creep-Fatigue 
Interaction 


Damage Mechanics Applied to 7/22 
Crack Initiation and Growth 
in Structures 


Large Space Structures: 8/27 

The Big Picture. Presenta- 
tion at 1st Workshop of the 
Mcgamechanics Research Con- 
sortium 

Static and Buckling Ana^y- 9/10 
sis of Composites 


Microstructure and Fracture 9/17 
Properties of Short Fiber 
Reinforced Thermoplastics 


Speaker (s) 

Leslie M. Lackman 
Director, Structures 
Rockwell International 
North TUnerican Aircraft Op. 

Robert G. Loewy 
Institute Professor 
Rensselaer Polytechnic Inst. 

Makoto Satoh 

Doctoral thesis dissertation 
Department of Mechanical 
Engineering, Aeronautical 
Engineering and Mechanics 

Jean LeMaitre 
Professor, University of 
Paris 6 , Laboratoire de 
Mdcanique et Technologie 
Cachan, France 

Robert G . Loewy 
Institute Professor 
Rensselaer Polytechnic Inst. 


N. G. R. Iyengar 
Department of Aeronautical 
Engineering 

I.I.T. 

Kanpur, India 

Dr. K. Friedrich 
Ruhr University 
Institute for Materials 
Bochum, West Germany 


Date 
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T ABLE VII>‘4 

COMPOSITES-RELATED V13ITS TO RELEVANT ORGANIZATIONS 
by RPI Faculty /Staff/Studonts 
(April 30, 1982 through SoptonUaer 30, 1982) 


Visited 


National Materials 
Advisory Board of 
the NAS , Washing- 
ton , DC 


Lewis Research 
Center, NASA 
(Workshop on 
Structural Dy- 
namics) 

SIGGRAPH '82 
Boston, MA 
(Panel on CAD/CAM 
Education) 

Vought Corpora- 
tion, Dallas/ 

Fort Wortli, TX 


Acoustic Emission 
Working Group, 
Knoxville, TN 


Date By 


5/26 Prof. S. S. 

Stornstein 


6/15 Prof. R. G. 
Loewy 


7/29 Prof. M. S. 
Shephard 


8/11, Prof. R. J. 
12,13; Diefendorf 
and 
9/2,3 

9/5-8 Prof. H. A. 
Scarton 


Purpose 

Participated in a 
special meeting on 
"Room Temperature 
Curing of Epoxy 
Resins" 

Presented talk; 
"Structural Dynamics 
Studies for Pan Tur- 
bine Designs" 


Presented talk: 
"CAD/CAM Education 
in the Engineering 
Curriculum" 

Discussed carbon/ 
carbon and ceramic 
matrix composites 
with B. Forscht 

Presented talk; 
"Acoustic Emission 
of Local Buckling in 
Thin-Walled Graphite 
Tubes " 
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coursQS £oir sovoral yoara. Bocauao of tho wide variety of 
special courses available throughout tho United States doal~ 
ing with this subject matter, only rather uniejuo aspects 
justify additional offerings in those areas. Our programs 
were planned from tho outset to bo unusual in tho respect 
that "hands-on" experiences wore inherent and required for 
completion of each course. Those aspects were enhanced for 
this third offering. Use of hand-hold programmable calcu- 
lators, application of optimissation programs using computer 
graphics, lay-up and cure of a simple part in graphite epoxy 
and testing tensile, shear and compression coupons to failure 
were all part of tho composite materials and structures 
course. As part of the second course, workshop sessions 
wore held in RPI's Computer Graphics Center. In those ses- 
sions tho enrolled engineers and managers defined, "con- 
structed" and analyzed their own finite element models on 
refresh graphics terminals. 

It was gratifying to note that NASA Army, Navy and Air 
Force, both largo and small industrial firms and other uni- 
versities were roprosontod among those who took tho courses 
and that their post-course questionnaires were overwhelmingly 
favorable . 

Faculty-staff and attendees for tho courses wore os 


follows ; 
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Composite Materials and Structures 
July 12-16, 1982 


Program Director; 

R. G. Loewy 
Institute Professor 
Rensselaer Polytechnic Inst, 

Faculty: 

R. Judd Diefendorf 
Professor, Materials Eng. 
Rensselaer Polytechnic Inst. 

H. Gunther Helwig 
Design Specialist 
Composites Group 
Dornier System Gmb H 

S . Leigh Phoenix 
Associate Professor 
Mechanical and Aeronautical 
Engineering 

Cornell University 

Students : 

Raymond Cannon 

Structural Integrity Engr. 

AFWAL/FIBE 

Wright Patterson AFB 
OH 45433 

Joseph A. Corrado 
Structural Engr. 

DTNSRDC 

Washington, DC 20590 

George Delli Santi 
Research Engr. 

Zimmer 

P. 0. Box 708 
V7arsaw, IN 46580 

Allan Ehlers 
Senior Staff Engr. 

The Bendix Corporation 
Guidance Systems Division 
Teterboro, NJ 07608 


Stephen W. Tsai 
Chief, Mechanics and Sur- 
face Interactions Branch 
Air Force Materials Lab. 
Wright-Patterson AFB 

Dick J. Wilkins 
Sr. Engineering Specialist 
General Dynamics 
Fort Worth 


Craig Golisano 
Mechanical Engr, 

Eastman Kodak 
901 Elmgrove Road 
Rochester, NY 14650 

Carl Holdren 
Design Engr. 

General Electric Company 
Interstate 75 and Neumann Way 
Cincinnati, OH 45215 

Karen Jackson 
Engineer 

U. S. Army Structures Lab. 
Langley Research Center MS 495 
Hampton, VA 23665 

A. Hugo Kruesi 
President 
Elfin Corporation 
550 Chippenhook Road 
Wallingford, VT 05773 
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Eugenio Rieghi 
Engineer 
S iai-Marche tti 
21018 Sesto Calende 
Via Indipendenda , 2 Italy 

Giantaolo Rusconi 

Engineer 

Siai-Marchetti 

21018 Sesto Calende 

Via Indipendenda, 2 Italy 

Dennis R. Schneider 
Aeronautical Engr. 
SA-ALC/MMSRE 
Kelly AFB 
TX 78241 


Thomas F . Smith 
Mechanical Engr. 
SA-ALC/MMETM 
Kelly 7VF3 
TX 78241 

Ray F. VJojcieszak 
Manager of Special Projects 
General Electric Company 
MD - J134 

Evandale, OH 45215 


The Finite Element Method 
August 2-6, 1982 


Program Director: 

M. S. Shephard 
Associate Director, Center 
for Interactive Computer 
Graphics and Assistant 
Professor of Civil Eng, 
Rensselaer Polytechnic Inst. 


Charles Ritter 
Vice President 
Jordan , Apostal , Ritter 
Associates Inc. 
Davisville , RI 

John Gray Harry Schaeffer 

Manager of Technical Services Pounder & President 
General Electric CAE, Inti. Shaeffer Analysis, Inc. 


Faculty; 

Richard H. Gallagher 
Dean, College of Engineering 
University of Arizona 


Students : 

Alfredo Cires 
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